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SECTION  1 
INTRODUCTION 

Hydrodynamics  has  in  recent  years  given  impressive  results  in 
the  solution  of  flow  field  problems  by  making  simplified  assumptions  of 
steady  flow  relative  to  a moving  body.  This  aspect  has  appeared  in  many 
physics  codes  attempting  to  solve  the  problem  of  the  rising  fireball. 
However,  because  of  computing  considerations,  flow  fields  have  generally 
been  ignored  in  large  system  codes  such  as  RANG  and  others.  This  has 
meant  errors  in  energy  deposition  and  a resulting  ghost  fireball  problem. 

A solution  to  these  problems  for  system  codes  is  presented  here. 

By  modeling  the  hydrodynamic  flow  using  the  simplified  assumption  of  steady 
flow,  the  detailed  particle  motion  can  be  specified  in  the  flow  field. 

This  model  is  a result  of  advances  that  have  been  made  in  the 
understanding  of  fluid  particle  movement  by  studying  the  deformation  of 
"drift"  of  material  surfaces.  Sir  Charles  Darwin  (1953)  has  drawn  impor- 
tant conclusions  about  the  movements  of  individual  particles  in  irrotational 
flow  of  fluid  around  hard  objects.  Darwin  considered  an  infinitely  thin 
plane  of  fluid  at  right  angles  to  the  motion  of  the  sphere  and  asked  what 
the  final  displacement  of  the  plane  was  after  the  passage  of  the  sphere. 

This  final  displacement  is  the  total  drift  function. 

Darwin's  method  has  been  used  to  calculate  the  position-time 
history  of  the  particles  of  fluid  as  they  move  around  a rising  fireball. 
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To  solve  the  rising  fireball  problem,  it  is  assumed  tliat  the 
fluid  is  incompressible  and  the  flow  is  irrotational , that  the  fireball  at 
any  instant  is  a perfect  sphere  with  no  entrainment.  This  reduces  the 
problem  of  finding  the  velocities  anywhere  in  the  field  outside  the  fire- 
ball to  one  of  taking  the  gradient  of  a velocity-potential  for  an  incom- 
pressible fluid,  the  motion  is  steady;  thus  the  fluid  does  not  cross 
steamlines. 

The  velocity-potential  (p  can  be  found  from  Laplace's  equation 
with  the  boundary  condition  of  no  radial  velocities  at  the  surface  of  the 
sphere  and  zero  velocity  at  infinity  for  a fluid  moving  around  the  sphere. 
The  velocities  and  positions  are  found  for  a coordinate  system  fixed  with 
the  sphere.  A Galilean  transformation  then  gives  the  position  in  a frame 
fixed  with  the  fluid. 


SECTION  2 

MATHEMATICAL  METHOD 

FLOW  FIELD  CALCULATION 

Consider  an  infinite  thin  plane  of  fluid  at  right  angles  to  the 
motion  of  the  fireball  vortex.  It  can  be  asked,  what  is  the  final  shape  of 
the  marked  fluid  plane  after  the  fireball  has  passed  through  it.  It  is  to 
be  expected  that  the  part  of  the  plane  nearest  the  spheroidal  fireball  is 
moved  the  greatest  distance.  The  fluid  contained  between  the  initial  plane 
position  and  its  final  position  equals  the  "hydrodynamic  mass"  or  "virtual" 
mass  associated  with  the  body's  motion. 

To  solve  the  problem  of  what  happens  to  the  displaced  fluid 
particles  requires  that  the  streamlines  be  determined  and  also  the  time 
at  which  each  fluid  particle  reaches  a given  point  measured  from  some  fixed 
reference  time  for  the  particles — in  this  case  when  it  passes  the  plane  at 
right  angles  to  the  center  of  the  spheroid  had  it  been  an  undisturbed  region. 

In  a coordinate  system  fixed  with  the  fireball,  with  z parallel 
to  the  direction  of  the  flow  and  the  reference  point  defined  as  2 = 0,  then 
we  require  solution  of  the  equations 

dt.  & (1) 

V V V 


^ M 
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t - z/u  -►0  as  z -*•  - 00 

where  is  the  z velocity  and  an  undisturbed  stream  flow  has  the  value 


V = u,  V = V = 0. 
z X y 


The  velocities  v , v , and  v are  found  from 
x’  y’  z 

V =’-^  V = V = 

X 3x  ’ y 3y  z 3z 


where  (}>  is  the  velocity-potential  solution  of  Laplace's  equation. 

This  definition  implies  that  far  upstream  material  planes  at  right 
angles  to  the  stream  are  planes  of  t=  constant.  This  is  the  classical 
"drift  function"  discussed  by  Darwin  (1953) . 

Consider  flow  past  a spherical  fireball  vortex  in  which  the  velocity 
field  far  upstream  from  the  obstacle  is  defined  as 


= V(x,y)  Vy  = v^  = 0. 


the  form 


The  stream  lines  of  the  flow  can  be  represented  by  equations  of 


X = x(xo,  yo,  z)  y = y(xo,  yo . z) 


where 


xo  = lim  (x) 
2 - 00 


yo  = lim  (y) 


2 - 00 


These  are  solutions  of  Equation  1.  The  solution  for  the  variable  t is 


t * t(xo,  yo,  z) 


1 ^ r ) ^ 

u J j V^(xo,  yo,  t)  u J 


dz 


(5) 


where  is  the  z component  of  the  velocity  on  the  stream  line  given  by 

Equation  3,  and  u is  the  undisturbed  fluid  velocity. 


Given  a point  in  the  flow  field  outside  of  a sphere,  the  drift 
function  can  now  be  determined. 


The  drift  function  t'  at  a burst  time  T'  can  be  found  by  sub- 
tracting the  actual  time  difference  AT  between  the  calculation  time  T and 
the  burst  time  T . Once  the  new  drift  function  t'  is  known,  a new  position 
can  be  found  from  it. 

Using  polar  spherical  coordinates  r,  6,  i)>  fixed  ac  the  center  of 
the  sphere,  the  Stokes  stream  function,  LAMB  (1932),  is  given  for  unit 
upward  velocity  as 

Pg  = r^sin^e  (1-aVr^)  (5a) 


where  a is  the  radius  of  the  sphere,  ij;  is  equal  to  zero.  It  should  be 
noted  that  po  x as  x -►  ®. 


The  solution  for  the  drift  function  t can  be  obtained  from  either 
of  the  coupled  ordinary  differential  equations 


Although  Equations  (6)  and  (7)  are  hyper-elliptic,  useful  expressions 
can  be  obtained  for  large  and  small  values  of  Po/a. 


For  large  Po/a,  Stokes  stream  function  can  be  expanded  in  powers 


of  a*/r^  as 


r = 


Po  /,  1 a^  • 3q  3 a*  . 6q  \ 

snfe(‘*  J -p  51"  9 - j ^ sin  6 . 


(9) 


Using  this  expansion.  Equation  (7)  can  be  expanded  to  give 


j.  Pod0  /,  3a®  . 6q  a*  . 9.  315  a*^  . i2q  \ 


(10) 


Integrating,  yields 


ut  = Po  Cot0  + 


|^Jsi„‘de-^Jsi„’dei|i|  ^Jsin-OdO 


..  (11) 


0 


0 


0 0 


The  limits  for  the  integral  satisfy  Equation  (8).  This  series  evaluation 
converges  for  all  0 with  the  value  po/a>  1.375. 


The  total  drift  function  z(P«)  from  Darwin  (1953)  is  given  as; 

:)  (12) 


z-*-( 


3 

al 

t^\ 

_ ^ 

315 

a”  j 

8 

Po 

\ 8 / 

P? 

'35/ 

128 

and 


It  should  be  noted  that  for  any  streamline 
t(0)  =2t(Tr/2)  -t(Tr-9) 


Z(po)  = lim  (ut-z)  = 2(ut) - lim  (ut-z) 
e-O  (0=7r/2  0^^ 


^^“*^^0=  Tr/2 


(13) 


For  small  values  of  Po/a  each  streamline  is  divided  into  two 
parts.  Part  one  for  0 near  0 or  tt  , and  second  one  on  which  (r/a  - 1) 
is  small. 

When  0 is  near  tt,  then  -sec0  * sin^0  in  Equation  (6)  and 
using  Equation  (5)  gives  the  value  of 


This  then  gives  the  result 


ut  = j Z(Po)  + 


Hn 


tan(y  6) 


2Po 

9a 


cot6  cosecd 


(16) 


The  value  of  the  total  drift  function  is  then  found  from  Equation  (16)  and 
(14)  to  be  for  small  (r/a  - 1)  to  be 


z(Po)  = 


5sfj)^ 


(17) 


Equations  (12),  (13),  (14),  (16),  and  (17)  can  be  used  to  calculate 
the  drift  function  t for  any  given  value  of  Po  , r and  9 . For  9 > S/bfr 
Equation  (14)  is  used.  For  1/2  tt  < 6 < 5/6  ti  and  r/a<1.5.  Equation  (16) 
is  a good  approximation.  When  r/a >1.75  Equation  (11)  is  valid.  The  gap 
between  1.5  < r/a<1.75  can  be  filled  by  interpolating  Equations  (16)  and 
(11). 


Lines  of  t = constant  drift  function  values  are  plotted  in  Figure  1. 
The  horizontal  lines  were  initially  at  right  angles  to  the  motion  of  the 
sphere.  They  are  distorted  around  the  sphere  as  it  moves  through  the  fluid. 
The  plot  was  made  in  a coordinate  system  fixed  with  the  rising  sphere. 

The  lines  of  constant  drift  function  show  in  detail  where  every 
particle  is  and  when  it  is  there.  For  example,  if  in  polar  spherical  co- 
ordinates, a point  is  given  as  r=1.6  and  9 = - tt/4  , the  corresponding 
drift  function  can  be  found,  which  is  say  + 2.0,  scaled  to  a unit  fireball 
with  unit  rising  velocity.  It  is  requested  to  know  where  the  point  was 
3.5  seconds  before.  This  would  put  it  on  a drift  function  plane  of  -1.5 
seconds,  with  the  same  stream  function. 
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Figure  1.  Plot  of  drift  functions  t and  stream  functions  ip  for  a 
unit  sphere  in  a unit  flow  field.  Coordinates  of  the  inter 
section  of  t and  \p  are  stored  in  the  table  XX  in 
Subroutine  CIPHER,  where  = W2  po^. 
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Jb- 


The  corresponding  new  position  at  a time  of  -1,5  seconds  can  be 
found  by  inverting  Equations  (11)  to  (17),  or  more  conveniently,  a table 
of  solutions  of  Equations  (6)  and  (7)  can  be  made  for  different  stream 
functions  and  drift  functions.  Since  both  the  stream  function  and  drift 
function  values  are  known,  the  solution  can  be  found  by  interpolating  from 
the  table.  Only  one  of  the  solution-variables  r or  9 need  be  stored 
since  the  other  can  be  calculated  from  the  constant  Stokes  stream  function 
Equation  (5a) . 

To  generate  the  table  that  inverts  Equations  (11)  to  (17),  the 
coupled  Equations  (6)  and  (7)  are  solved  numerically  in  the  form: 


d0  _ /,  . a \ sin0 

dt  " ■ \ 2r^l  r 


where  the  rise  velocity  and  radius  a are  both  taken  as  unity. 


(18a) 

(18b) 


Numerical  integration  of  Equations  (18a)  and  (18b)  gives  the 
position  of  the  particle  in  the  flow  field  for  discrete  times.  The  numerical 
method  used  was  that  of  Gear  (1971).  The  Gear  algorithm  is  useful  for  solving 
stiff  equations.  It  uses  an  Adams  Predictor-Corrector  method  with  the  initial 
values  generated  by  a Runge-Kutta  method. 


The  polar  spherical  radial  coordinate  r was  saved  in  a data  table 
for  each  discrete  time  t and  stream  function  \p  , where 


* j r^sin^0 


(19) 
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Thus  for  a given  ^ and  t,  r can  be  found  by  a table  look-up. 
Then  9 can  be  found  from  the  inversion  of  Equation  (19)  using  r and  ip 


0 = sin  * 1 -4- 


Thus,  the  position  is  specified.  The  table  look-up  and  interpolation  is 
done  in  subroutine  CIPHER  and  the  outputs  converted  to  fireball  centered 


Cartesian  coordinates. 


Once  the  positions  (z^,  x^)  in  fireball  coordinates  is  known,  a 


simple  Galilean  transformation  gives  the  position  (Z^,  X^,  Y^)  in  a coordin- 


ate system  fixed  at  the  position  of  the  initial  burst. 


Zc  = Sri  - V + z 

f R q 


where  Sri  is  the  slant  range  between  the  initial  burst  position  and  the 


present  fireball  position,  V is  the  fireball  rise  velocity,  and  z is  the 

K q 


z coordinate  fixed  at  the  initial  burst  position.  Coordinates  and 


remain  constant  under  the  Galilean  transformation  in  the  z direction. 


The  point  defined  by  X^,  Y^,  and  , in  the  initial  burst  position 
coordinates,  is  then  rotated  to  the  earth  centered  Cartesian  coordinates. 
The  vector  position  in  this  new  frame  is  then  added  to  the  vector  position 
of  the  initial  burst  giving  the  new  vector  position  of  the  point  in  earth 


centered  Cartesian  coordinates. 


The  trajectory  of  the  particles  in  a frame  fixed  with  the  stationary 
fluid  are  plotted  in  Figure  2 for  some  typical  values  of  Po(=  .5,  1.0,  1.5, 
1.8).  From  the  figure,  it  can  be  seen  that  as  the  fireball  passes  the  point, 
it  is  at  its  maximum  lateral  displacement.  It  should  also  be  noted  that  as 
the  fireball  passes  the  point,  the  point  takes  on  a retrograde  motion. 
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SHOCK  DISPLACEMENT  CALCULATIONS 


The  shock  displacement  is  calculated  after  the  point  has  been 
moved  back  in  time  to  a burst  time.  It  is  assumed  that  the  shock  arrives 
instantaneously  at  the  point  of  interest.  This  is  a zeroth  order  approxi- 
mation and  will  be  changed  to  include  the  shock  arrival  time  and  shock 
duration.  For  calculation  times  long  enough  after  burst  time  for  shock 
traversal  of  the  point,  this  assumption  causes  only  a small  error  because 
of  two  facts:  1)  if  the  fireball  is  very  far  away  from  the  point,  the 
shock  displacement  time  is  greatly  in  error,  but  the  actual  displacement 
is  very  small,  2)  if  the  burst  is  close  to  the  point,  the  timing  error 
is  small  and  the  shock  displacement  is  large. 


The  distance  from  the  initial  burst  position  to  the  point  is 
scaled  to  dimensionless  coordinates  using  a modified  Sach's  scaling 


p 1/3  Si/3H 

(''ktp7)  " 

where  is  the  yield  in  kilotons,  P 

interest,  Po  is  the  sea  level  pressure, 
above  (+)  or  below  (-)  the  fireball,  H 

s 

expression  reduces  down  to 


is  the  pressure  at  the  point  of 
Si  is  the  height  of  the  point 
is  the  scale  height.  This 


Si/3H 

a = 13.6  e ^ 


where  is  the  initial  fireball  radius  at  the  end  of  the  radiation 

phase. 


Based  on  a fit  to  a RADFLO  run  (5KT  at  9.15  km),  a shock  displace- 
ment AA  in  scaled  distance  is  given  by 

AA  = 0.096  (A'  (2+A"))‘‘ 

where  A"  = A + AA  is  the  final  position  in  scaled  units. 
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SECTION  III 


CODING  CONSIDERATION  FOR  FLOW  FIELD  CALCULATIONS 

CONTROL  OF  CALCULATIONS  BY  SUBPROGRAM  HYDRO 

The  above  mathematical  method  of  calculating  particle  movements 
in  the  flow  field  around  a rising  fireball  has  been  coded  into  eight  sub- 
programs for  use  in  ROSCOE.  The  controlling  routine  is  HYDRO,  which  is 
called  with  a given  time  and  position. 

Previous  to  a call  to  HYDRO,  a call  has  been  made  to  the  burst 
environment  module  (PMDS)  which  sets  up  the  vortex  radii  at  the  initial 
times  and  at  the  calling  time.  Also  saved  are  all  the  fireball  vortex 
positions  at  that  time  of  call. 

Given  the  position  of  a point,  the  time  and  pertinent  fireball 
positions,  HYDRO  makes  a call  to  subprogram  EDITX  to  determine  which  bursts 
affect  the  given  point. 

EDITX  begins  the  loop  over  all  existing  fireballs  including  those 
created  by  a hydromerge.  If  the  burst  was  created  by  a hydromerge,  it  creates 
the  initial  vortex  information  for  the  merged  fireball.  It  then  calculates 
the  distance  of  closest  approach  to  the  line  connecting  the  initial  fireball 
center  and  present  center.  If  this  radius  scaled  distance  is  greater  than 
1.8,  it  rejects  this  fireball  for  consideration.  If  the  time  it  takes  for 
the  fireball  to  get  to  the  point  is  very  large,  the  event  is  rejected.  If 
the  point  is  more  than  a radius  on  front  of  the  fireball,  the  burst  is  re- 
jected. For  all  bursts  accepted  as  influencing  the  given  point,  a count  is 
kept  (=NAI)  and  the  fireball  index  is  stored  in  the  NCAG  array. 
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HYDRO  then  begins  to  loop  over  all  the  bursts  including  bursts 
created  during  a hydromerge.  It  starts  with  the  most  recent  burst  time 
first  and  does  the  calculation  back  to  this  time.  It  then  does  the  shock 
displacement  calculation.  It  saves  this  new  calculated  position  and  puts 
it  into  the  output  array  if  the  time  was  a burst  time  and  not  merely  a- 
merge  time.  The  saved  new  position  is  now  used  to  regress  the  calculation 
back  to  the  second  to  the  last  occurring  burst.  This  new  position  is  saved 
i and  used  to  continue  the  calculation  back  to  where  the  parcel  was  for  the 

[ first  burst  before  the  arrival  of  any  shock  waves. 

, For  each  regression  back  to  a burst  time,  an  inner  loop  is  made 

over  all  the  fireballs  that  will  influence  the  point  during  this  time 
increment.  After  the  time  regression  goes  back  beyond  the  burst  time  of 
an  event  affecting  the  point,  this  event  is  removed  from  the  NCAG  array, 
and  the  calculation  continues  with  the  remaining  influencing  fireballs. 

When  multiple  events  influence  the  point,  the  individual  dis- 
placements from  each  burst  are  saved.  After  the  inner  loop  over  all 
influencing  events  is  completed,  a vector  sum  is  made  of  all  the  individual 
displacements. 

A check  is  made  of  the  influencing  fireball  type  flag  KINDF  to 
determine  if  the  fireball  is  active  at  this  time.  If  it  is  not  active, 
i.e.,  has  been  radiation  or  hydromerged,  then  it  is  not  used  to  calculate 
the  change  in  position  of  the  point  PXYZ.  After  the  calculation  regresses 
back  to  the  time  of  the  merge  for  this  fireball,  the  burst  is  used  to 
determine  the  position. 

‘ i 
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CRITERIA  FOR  EDITING  FIREBALLS 

Before  a flow  field  calculation  is  performed,  a call  is  made  to 
an  editing  routine.  This  routine  loops  over  all  the  existing  fireballs, 
active  and  inactive,  and  determines  which  fireballs  will  affect  the  given 
point.  All  fireballs  that  will  ever  have  an  affect  on  the  point  are  stored 
in  array  NCAG. 

If  the  distance  of  closest  approach  to  the  line  of  centers  between 
the  initial  burst  and  present  fireball  position  (scaled  to  the  fireball) 
radius  at  the  calculation  time)  is  greater  than  1.8,  the  fireball  is  rejected 
as  not  affecting  the  point. 

If  the  point  of  interest  is  more  than  1.5  radii  above  the  fireball 
position  at  calculation  time  or  1.5  radii  below  the  initial  burst  position, 
the  burst  is  rejected. 

SINGLE  BURST  CASE 

For  the  single  burst  case,  when  EDITX  returns  NAI  equal  to  1 and 
NCAG(l)  equal  to  1,  subprogram  HYDRO  begins  the  calculation  of  the  effect 
of  burst  1 on  PXYZ.  The  burst  will  always  be  active  since  there  is  only 
one  burst.  Also  the  calculation  can  proceed  back  to  the  initial  burst  time 
in  one  pass  through  the  routine. 

Hydro  calculates  the  distance  of  closest  approach,  DIST.  It  cal- 
culates the  slant  range  along  the  line  of  centers  between  the  initial 

burst  position  and  the  position  at  calling  time,  the  slant  range,  SI,  from 
burst  point  to  point  of  closest  approach;  slant  range,  from  present 

burst  position  to  point  IT  the  point  is  more  than  a fireball  radius 

below  the  initial  burst  point,  the  effect  of  burst  N on  the  point  is 
considered  negligible. 
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The  radius  of  the  fireball  where  most  of  the  influence  is  occurring 
on  point  P is  found  by  interpolating  linearly  between  the  initial  and 
final  radii: 

R = (R.  S . + R S-)/S  . 

^ in  rd  V V rl 

where  is  the  slant  range  from  the  present  position  of  the  fireball  to 

the  point  of  closest  approach,  R^^  is  the  initial  vortex  radius,  R^  is 
the  present  vortex  radius. 

MULTIPLE  BURST  CASE 

The  mathematical  method  for  the  case  where  two  or  more  fireballs 
are  affecting  a point  is  handled  similarly  to  the  single  burst  case.  The 
displacement  calculation  is  done  separately  for  each  fireball.  The  final 
displacement  is  the  vector  sum  displacement  over  all  fireballs. 

No  attempt  is  made  to  account  for  the  change  in  a fireball  flow 
field  due  to  another  fireball  flow  field.  It  is  expected  that  the  effect 
is  small  except  in  cases  where  the  fireballs  are  closer  than  .5  fireball 
radii. 


SECTION  IV 

DESCRIPTION  OF  SUBPROGRAMS 

Subroutine  HYDRO  is  the  main  ROSCOE  subprogram  for  calculating 
the  particle  time  history  in  the  flow  field  around  a rising  fireball. 


Fireballs  are  edited  as  to  their  effect  on  the  given  point  by 

EDITX. 


Subprogram  SYZYGY  calculates  the  drift  function  t (=TOUT)  for 
the  initial  given  point. 

Subprogram  CIPHER  calculates  the  position  of  the  particle  at  a 
given  drift  function  and  stream  function. 

Subprogram  SCHCK  calculates  the  shock  displacement  of  a given 
point  by  a given  burst. 
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Description  of  Code  HYDRO 


Purpose 

Subroutine  HYDRO  calculates  the  position  of  a given  point  at  a 
given  time  at  all  burst  times  previous  to  the  given  time.  This  requires 
that  the  code  run  backwards  in  time. 

Inputs  to  Subroutine  HYDRO 

Formal  Argument 


PXYZ(I)  = 3 word  array  of  coordinates  x,  y,  z in  cm. 

Earth  centered  Cartesian  with  x through  Greenwich 

Inputs  from  Labeled  Common 


Block 

Variable 

Description 

/RADDAT/ 

XIN(3,I) 

Initial  vector  position  of  event  I (cm) 

XEV(3,I) 

Vector  position  of  event  I at  time  = TSIM  (cm)  j 

REIN(I) 

Initial  vortex  radius  =4*  RHZERO  j 

REV(I) 

Vortex  radius  of  event  1 at  time  = TSIM  (cm)  I 

VRZ(I) 

Event  rise  velocity  magnitude  (cm/sec)  1 

/EVXDAT/ 

NUMX 

Number  of  real  events  ' 

TSIM 

Time  of  calculation  (sec)  i 

/EVENTX/ 

NF 

Total  number  of  events  at  time  = TSIM 

< 

includes  hydromerges  j 

TB(I) 

Burst  time  of  event  I 

/GEOTD/ 

KINDF(I) 

Flag  to  indicate  type  of  event  j 

= 1 spheroidal 

= 2 spheroidal,  pressure  equilibrium 
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3 torus 


= 4 merged  during  radiation  phase 
= 5 merged  during  hydrodynamic  phase 
MRGID(I)  Flag  to  indicate  index  of  merged  events 

TCHAR(I)  Characteristic  time  for  merged  events 

/SAVEVX/  BUFFR(15,I)  Initial  event  radius  for  the  Ith  event  to 

which  all  parameters  arc  scaled  (=RHZERO(I)) 
in  common  block/EVENTX/ 

Output  from  Subroutine  HYDRO 
Formal  Arguments 

XARAY(3,N)  Vector  position  array  for  point  PXYZ  at  all  previous  events 
(N=1,NVMX)  (cm) 

MODE  Flag  to  indicate  significant  fluid  movement 

= 1 fluid  has  moved 
= 0 fluid  has  not  moved  significantly 

External  Subprograms  Used  in  HYDRO 
Name  Information  requested/description 

EDITX  Loops  over  all  bursts  and  selects  those  that  will  affect 

point  PXYZ 

SYZYGY  Gives  the  time  since  the  point  passed  the  center  of  the 

burst  (+)  or  the  time  it  will  take  to  pass  the  center  of 
of  the  burst  (-),  scaled  to  radius  of  fireball 

CIPHER  Calculates  the  position  of  the  point  for  a given  time  in  a 

coordinate  system  attached  to  the  burst  center  and  in  the 
x-z  plane.  It  does  the  calculation  for  a unit  vortex 
radius  and  unit  velocity. 

EVENAD  Performs  the  vector  addition  of  point  displacement  resulting 

from  the  multiple  burst  interactions. 


SCHCK  Calculates  the  shock  displacement  from  each  burst  on  the 

point  at  the  burst  time, 

ROOTT  Calculates  the  time  at  which  a given  stream  function  crosses 

the  zero  axis. 

DISCAP  Calculates  the  distance  of  closest  approach  of  the  point 

to  the  line  connecting  centers.  Calculates  the  needed 
slant  ranges. 

External  Utility  Subprograms  Required 


Name 

Vector  package 

LOCLAX 

Generates  transformation  matrix 

XMIT 

Copies  vector  A to  B 

XMAG 

Gives  magnitude  of  vector 

SUBVEC 

Subtracts  two  vectors 

VECLIN 

Adds  two  vectors 

MATMULT 

Does  matrix  multiple 

DOT 

Does  Dot  product 

UNITV 

Finds  unit  vector 

CROSS 

Finds  cross  product  A^B 

CROSl 

Find  cross  product  and  normalizes  to  1 

PROJ 

Finds  projection  of  A on  B 

LOCI 

Find  SCM  address 

FDIV 

Prevents  division  error  abort 

VECSUM 

Vector  sum 

EFCGEO 

Transforms  from  xyz  to  geographic  coordinant 

GEOEFC 

Inverse  of  EFCGEO 

L . . _ J 


GLOSSARY  OF  VARIABLES 


NAME 

PXYZ(3) 

XARAY(3,I) 

MODE 

TIM 

TSIM 

XSAV(3) 

Ml 

L 

NAI 

MiA(I) 

NRUN 

NX 

K 

XOU(3) 

KFF 

TOBIO 

TBCK) 

DTEL 

NCAG(IO) 

NN 

KINDF(IO) 
TCHAR(IO) 
PS  I 
VNORM 


Roscoe  Routine  HYDRO 


Description 

Input  vector  position 

Output  vector  position  at  each  burst  time  of  Index  I 

Flag  to  indicate  significant  flow  field  movement 

Calculation  time  local  variable 

Global  calculation  time 

Local  Var.,  vector  position  array 

Flag  to  indicate  data  is  saved  for  fireball  vortex 

Running  index  for  burst  index  contained  in  NCAG 

Return  from  EDITX,  Number  of  burst  affecting  PXY2 

Array  of  flags  indicating  data  saved  for  burst  I 

Dummy  running  variable  for  DO  1000 

Number  of  bursts  up  to  time  TIM 

Index  of  burst  starting  with  last  burst  first 

Local  variable,  array  of  vector  position 

Flag  whether  index  K is  for  real  burst  or  merged  burst 

Local  variable  set  to  the  burst  time  or  merge  time 

Burst  time  from/GEOTD/ 

Time  difference  between  present  time  and  burst  time 
Array  of  burst  indexes  affecting  PXYZ 
Running  index  of  burst  affecting  point  PXYZ 
Array  of  flags  describing  fireball  type/GEOTD/ 
Characteristic  burst  or  merge  time  for  each  fireball 
Stream  function 

Speed  of  fireball  in  effective  vortex  radii/sec 
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TZ 

SRI 

SI 

DIST 

R 

XIN(3,I) 

XEV(3,I) 

RAD 

REV(I) 

SRD 
RH 
RR 
PS  I 

VRZ(I) 

TRELX 

XF 

YF 

ZF 

XQ 

ZQ 

D(3) 

CVEC(3,N) 

UVEC 

NRA 

NCARB 

NCB 

P0F(3) 

NUMX 


Relative  time  since  passage  of  point  past  fireball  (normalized 
by  R) 

Slant  range  of  initial  position  to  present  position 

Slant  range  from  initial  position  to  PXYZ  on  line  connecting 

centers 

Distance  of  closest  approach  of  PXYZ  to  line  of  centers 

Radius  of  vortex  near  PXYZ  used  to  scale  distances 

Initial  position  of  fireball  vortex 

Final  position  of  fireball  vortex 

Slant  range  from  fireball  position  to  PXYZ 

Radius  of  vortex  at  calculation  time 

Slant  range  from  fireball  to  point  of  closest  approach 

Normalized  distance  of  closest  approach 

Radius  normalized  RAD 

Stream  function  in  unit  velocity  field 

Speed  of  fireball  unnormalized  (cm/sec) 

Time  difference  scaled  to  fireball  with  unit  rise 

X position  in  fireball  centered  coordinates  of  point  PXYZ 

Y position  in  fireball  centered  coordinates  of  point  PXYZ 

Z position  in  fireball  centered  coordinates  of  point  PXYZ 

X positions  of  point  at  new  time  TIM 
Z positions  of  point  at  new  time  TIM 

Vector  position  of  point  PXYZ  in  initial  fireball  coordinates 
multiply  used 

Array  of  vectors  from  initial  fireball  N to  final  fireball  N. 

Transformation  matrix  from  XYZ  to  user  frame 

Row  size  of  matrix 

Column  size  of  matrix 

Vector  column  size 

Vector  position  in  burst  centered  coordinates  at  time  = TOBIO 
Number  of  actual  bursts 
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SP0U(3,N) 

TMBU(IO) 

NCOT 

ALPHA 


Local  array  to  store  positions 
Burst  times  in  descending  order 
Total  number  of  print  times 
Shock  scaling  parameter 


I— J f!»IO  SlMT  RANQ($ 
I SI.  SHI  MO  OJST 


Flow  Diagram  for  Subprogram  HYDRO 


TMBU(NCOU) 

TOBIO 

EVENT  TIME 


SAVE  CALCULATED 
POSITION 
XOU  • SPOU 


This  subroutine  loops  over  all  existing  fireballs  at  calling  time 


and  determines  which  fireball  will  affect  the  given  point. 
Inputs 

PXYZ(3)  Vector  position  of  point 

Inputs  from  Common 

XIN(3,L)  Initial  position  of  fireball  L 
XEV[3,L)  Final  position  of  fireball  L 
TB(L)  Burst  time  for  event  L 


Outputs  to  Common 

NAI  Number  of  bursts  interacting  with  point 

NCAG(I)  Array  of  burst  indexes  interacting  1=1,  NAI 

Externals  used 

DISCAP  Calculates  distance  of  closest  approach 

Method 

EDITX  is  a simple  routine  which  is  called  every  time  HYDRO  is 

called. 

It  loops  over  all  fireballs  and  calculates  the  distance  of 
closest  approach  of  the  line  connecting  centers  with  the  given  point.  The 
calculation  is  done  in  DISCAP. 
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Criteria 

A point  is  rejected  if  it  is  greater  than  1.8  final  fireball  radii 
from  the  line  connecting  centers. 

It  is  rejected  if  it  is  1.5  fireball  radii  above  the  final  posi- 
tion or  1.5  radii  below  the  initial  position.  The  radius  at  call  time  is  used. 

Glossary  of  Variables 
Name 

PXYZ(3)  Input-vector  position 

NAI  Output-number  of  burst  affecting  point 

NCAG(L)  Output  - array  of  burst  indexes 

L Running  variable  over  all  bursts 

NX  Total  number  of  bursts 

SR12  Slant  range  from  final  fireball  position  to  initial 

CVEC(3,L)  Vector  from  initial  fireball  position  to  final  position 

XIN(3,L)  Initial  position  vector 

XEV(3,L)  Final  position  vector 

DIST  Distance  of  closest  approach  to  line  of  centers 

SI  Distance  from  initial  position  to  point  of  closest  approach 

TDIS  Time  for  fireball  to  travel  from  point  of  closest  approach  to 

intiial  position 

REV(L)  Final  vortex  radius  (event  L)  (cm) 

VRZ(L)  Speed  of  rise  for  event  L (cm/sec) 

TTEST  Time  to  go  from  point  to  1.5  fireball  radii  below  initial 

fireball  position 
NCO  Running  variable 
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Description  of  Code  SYZYGY 
Purpose 


i 
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I 
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This  subroutine  calculates  the  drift  function  TZ  scaled  to  a unit 
fireball  radius  rising  at  unit  speed,  given  a position. 

Inputs 

R Radius  of  fireball 

VR  Rise  speed 

PSI  Stream  function  (=  1/2  Po) 

Outputs 

TOUT  Drift  function 

Outputs  to  Common 


THETA 

Angle  to  point  measured  from  rise 

line 

RR 

Distance  in  radii  R from  fireball 

to  point 

RH 

Distance  of  closest  approach  in  R 

VNORM 

Speed  scaled  by  R 

XROZ 

Total  drift  displacement 

TZ 

Drift  function  time  (sec) 

Method 

See  Section  II 


'1 
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Glossary  of  Variables 
Name 

COSEC  Cosec  of  THETA 

t 

i COTAN  Cot an  of  THETA 

[ 

ROA  Equivalent  to  RO 

RO  = Po 

XROZ  Total  drift  function 

UTZ  Distance  from  reference  (z=0) 

All  other  variables  are  defined  previously. 
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Description  of  Code  CIPHER 
Purpose 

This  subroutine  calculates  the  position  of  a point  in  fireball 
coordinates  for  a given  stream  function  and  drift  function. 

Inputs 

T 
PP 

Inputs  from  Common 
THETA 
VNORM 
RH 
RR 

Outputs 

X 
Y 

Method 

The  coupled  equations  (6)  and  (7)  in  Section  II  were  solved  and  the 
r coordinate  solution  was  stored  in  a table  for  35  times  and  15  stream 
functions. 


Drift  function 
Stream  function 

Angle  to  point  measured  from  line  of  rise 
Rise  speed  in  radii 
Distance  of  closest  approach  in  R 
Slant  range  from  fireball  to  point 

Cartesian  coordinate  perpendicular  to  rise  in  R 
Cartesian  coordinate  parallel  to  rise  in  R 


Then  for  a given  time  (drift  function)  T and  stream  function 
PP,  a value  of  r can  be  picked  from  the  table. 
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The  X and  Y coordinates  are  now  determined  from  the  equation 


and 


/2*PP*r^ 

J r^-r 


y = 


The  sign  of  Y is  positive  above  the  fireball  center  and  negative 
below.  Close  to  fireball  the  sign  of  Y is  determined  by  finding  the  slope 
of  the  stream  function  dy/dx.  For  values  of  X far  from  the  sphere 
function  subprogram  ROOTT  is  used  to  find  the  time  at  which  the  given 
stream  function  crosses  the  zero  Z axis. 
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Glossary  of  Variables 
Name 

T Drift  function 

PP  Stream  function  (=i/^Po) 

X Output  - Cartesian  coordinates  perpendicular 

Y Output  - Cartesian  coordinate  parallel  to  rise 

NFLAG  Flag  to  indicate  if  point  has  moved 

TI  (35)  Array  of  drift  functions 

PS(15)  Stream  functions  (=i/2p5) 

XX(35,15)  Array  of  solutions  for  TI  and  PS(=r) 

P Local  variable  - stream  function 

NOFF  Flag  indicates  whether  P is  off  table 

ITOFF  Flag  indicates  whether  T is  off  table 

FT  Interpolation  weight  for  T 

FP  Interpolation  weight  for  P 

XN  Nth  point  to  interpolate  from 

RQ  Interpolated  solution  r 

XTEM2  X^ 

YTEM2  r^-x^ 

INFLG  Flag  indicates  whether  this  is  first  or  second  pass  through 

calculation 

XS  Temporary  x coordinate 

DELX  Change  in  X from  pass  1 to  2 

DTEL  Time  difference  (sec) 

TDIF  Time  difference  since  it  made  it  to  edge  of  table 

DNOM  x^-2«P 


Description  of  Code  SCHCK 
Purpose 

This  subroutine  calculates  the  shock  displacement  of  a given 
point  from  a given  fireball 

Inputs 

X0U(3)  Initial  position  vector 

K Fireball  index 

ALPHA  Scaling  parameter 

Inputs  from  Common 

XIN(3,L)  Initial  position  for  fireball  L 
XEV(3,L)  Final  position  for  fireball  L 


Output 

XSAV 


Vector  position  after  displacement 


Method 

See  Section  II.  It  does  not  take  into  account  the  change  in 
the  fireball  position  resulting  from  the  interaction  of  one  fireball  with 
another. 


A more  sophisticated  shock  routine  should  take  into  account  the 
change  in  position  of  all  other  fireballs  from  the  shock  passage  of  the 
fireball  under  consideration.  This  however  would  slow  down  the  calculations 
considerably  by  necessitating  a NxN  calculation  of  shock  interactions  and  by 
causing  a re-evaluation  of  flow  field  parameters  for  every  burst. 


36 


1 


Glossary  of  Variables 
Name 

C(3)  Vector  from  fireball  to  point 

SR  Magnitude  of  C 

SRLAM  Scaled  slant  range 

DLAM  Scaled  shock  displacement 

SRNEW  Final  slant  range 

XSAV  Final  position  after  shock  passage. 

All  other  variables  are  defined  previously. 
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Varriables  Generated  by  MODEL  and  used  by  ROSCOE  model  HYDRO 

Common  Block  Array  Used 

Block  Variable  Description 

GEOTD  TCHAR(I)  Array  of  characteristic  times  for  event  I 

KINDF(I)  Type  of  event  I,  i.e., 

= 1 spheroidal 

= 2 spheroidal  pressured  equilibrium 
= 3 torus 

= 4 merged  during  radiation  phase 
= 5 merged  during  hydrodynamic  phase 
MRGID(I)  Flag  to  indicate  merging 

For  non-merged  events  CMRGID(I)=KINDF(I)) 

For  merged  events,  MRGID(I)  units  digit  contains 
the  event  index  from  which  the  merge  took  place; 
the  tens  digit  the  event  index  of  the  second  event. 
The  MRGID  of  the  events  from  which  the  merged  event 
occurred  contains  in  the  unit  digit,  the  other  event, 
in  the  tens  digit  the  new  event  formed. 

RADDAT  XIN(3,I)  Initial  geocentric  Cartesian  coordinates  with  X 

through  Greenwich  for  the  event  I,  with  ordering 
X,Y,Z,  (cm). 

XEV(3,I)  Geocentric  Cartesian  coordinates  of  event  I at  the 
calling  time 

REIN  (I)  Initial  vortex  radius  (cm) 

REV(I)  Vortex  radius  at  calculation  time(cm) 

SAVEVX  BUFFR(15,I)  Initial  event  radius  for  the  event  (cm)  to  which 

all  parameters  are  scaled  (=RHZER0(I))  in  common 
block/EVENTX/ 
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The  output  array  XARAY  is  filled  in  time  ascending  order. 


Running  times  for  multi-event  scenarios  at  relatively  high 
altitude  and  high  energies  promise  to  be  long.  A five  burst  scenario 
with  a hydromerge  and  radiation  merge  for  example  will  run  greater 
than  77  ms  per  call  with  multiple  events  affecting  the  point.  (Run 
on  a CDC  7600) . 
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SECTION  V 


SUBPROGRAM  LISTINGS 


myobo 


SUbROUTImE  HYn»P(P*V2,X»R»Y,MOOE) 


TMiS  SUHPROGb**'  C*LCUL*TE5  THE  Tint  hISTURY  OF  THf  GIVEN  POINT 
PXyZ  *nO  returns  ITS  position  at  each  EVENT  TiHEt 


INPUT 

pXyZSARRtV 


Of  GEOCENTRIC  CARTESIAN  COUHOINaTES- 


OUTPuT 

X*RAT(5»I) 


VECTOR  POSITION  OF  EVENT  I AT  EVENT  Tl"E, 


COMHON/EvXOAT/NUHX,NAl,XAR(3,10),NCAG(lO),TSiri,NCOT,THBU(10) 


HYDRO,? 
•hydro. S 
hydro. a 
hvdPO.S 

HYDRO, N 
hydro,; 
hydro, 8 
-hydro, 9 
hydro, lO 
hydro, It 
hydro, t? 
hydro, iJ 
hydro, I# 
hydro, j 8 
evyoat.? 
f v»OaT. 3 

EVENTX.? 

Evfntx.3 


EVENTX  PaRAMITERS 

COpHON  /EVENtx/  NX,  lOX,  TB(GO),  HtatbD),  CLBlbO),  GEBtbO), 

1 lOGAO(SO),  RmOHCjO),  hSHISO),  TE-BtbO),  V R I SE  ( 50  J , E Vt  N TX> 

-? RD2EOO(50),  HHZEROtSu),  RUZERO(So),  BXb(50), EVEnTX.S 

3 eVBIbO),  BZBtSO),  LHVb(SO) , XAlPhA (SO) , nALCM  EVfNTX.b 

fQ|)IVAEENCE(TCH(t)»TE)(l))  . hydro, I8 

TIhe  dependent  MnDEL  PARahETERS  CFDTD,? 

COmhON  /RFOTd/  nF,  InOXF(SO),  KTF(SO).  RLKSO),  hFCSO),  CCF(50),  r.FoTU.3 
1 ClFISO),  HhXXKSO),  hhinE(BO),  KINCKSO),  IlLTK50),&tf.TD,a 

-2 AGEtSO),  NoTA,  I\0*0(100),  UEAoL(tOO),  nOR(IOO),  CEOTU,? 

3 HoBdfto),  RTbstlon),  PL’tS(ino),  HVsdno),  NM)SdOO),&FoTC,h 


« GCBTAd0(l),CLKrA(10O),Tfi«;0),TEHAW(S0),MWCl0(S0), 

5 XFR(50)iYFH(SO),Z)fl(SO),RUT(50) 

OIhFNSiPn  TOBd) 

COMHON/RADDAT/XIN(i,lO),X£V(l,lO),RElNdO),HEV(10),VXYZ(3,lO), 

— I VRZ(io),eVEC<J»tO) 


COhHON/GeOMO/THETA,RR,RH,VNORh,DTEL,XHOX,TI. 


MATHehaTjcal  XNO  geophysical  constants  .. 

COmMON  /cnSTnT/  re,  pi,  hALFPI,  RaOIAn,  DEGPRD,  CRt? 
DIpEnSjCn  P«YZC3)»XSAVC3),X0Ut3)iUytL(3,3)»C{3),0(3)»P0f 131.- 
t Spr,|(i,inj,XAHAY(3,lO) 

COp*'ON/FiC'-/PSlA(tO),VNCRHAdO),TZAdO),ELYAdO),AZPTAdO), 

I SRlA(t0),81*(t0),DlSTAd0),RAd0),mAd0) 


c-. 

c 

c 

c 


0*T*  TRAVL/1.E05/ 


call  up  the  edit  routine 
FRoh  being  considered  in 


CtnTD.T 

GEOTD.S 

HYDRO, go 

RACDAT.? 

raddat.i 

RADOAT.a 

GECO,? 

CEOhd.j 

CNSTNT,? 

CnsTNT.3 

HYDRO. ?<! 

hydro, ?S 

hydro, ?6 

HYDRO,?? 

hydro, ?8 
HYDRO, ?9 
hydro, 3O 

HYDRO, J! 


6 

11 

13 

1« 

IS 

?« 

?s 

- ?6 


30 


TIh»T81“ 

C*lU  XhIt(5,pxYZ,X8Av)  ■ 

C*iL  t0lTX(X8AV) 

Mlat 

L*NAl 

DO  30  IBl.tO 
Hi  A ( I ) bO 

continue 

nCoUbO 


TO  ELIhINaTE  The  nOn  essential  lvents 
the  FLO-  PILLO  calculations  FOR  THE  POInThydRQ.I? 

HYDRO, 3 3 


MrDH0.3« 
hydro, 35 
. HYDRO. 3S 
hydro, 3? 
hydro, 38 
hydro, 3R 

HYDRO, aO 

hydro, «1 

HYDRO, u? 


1 


40 


HTOHO 


so 


52- 

— . 

35 

56 

99 

90 

96 

50 

50 

■—  51 

52 

60 

59 

55 

80 

55 

5T 

61 

- - 

62 

65 

lOo 

65 

c 

C 



C-  - 

67 

72 

75 

76 

105 

76 

100 

c 

c 

c 

c 

c 

102 



105 

109 

106 

tor 

til 

— 112 

.. 

11« 

116 

116 

Ho 

116 

ISO 


no  1000  NRUN«t>NX  HVORQ.aS 

xsmX"NBUN*l  — • HvnBC.ao 

LOoP  CV€r  fWENTS  INClUOINS  hYOROmERGLO  evfNTSi  ►"'ORO.dS 

hydro. u* 


begin 

00  MOST  hfCENT  event  first 

C*LL  Xh1t{3»*S*V*1«0U) 

KFr»o 

lFfKlNDF(K),GT.S)  CO  TO  60  — 
lFtHhGlD(F),efj,K)  Go  TO  60 

TOBlOaTCH4R(K)  

kFfsI 

GO  TO  80 

TObI08T8(k) 

kFfsO  

continue 

(iTfL  *TIM-T09i0  

iFfSAI.LE.O)  RO  TO  TOO 

iTrNCt&m.Eo^O)  GO  TO  700 

iF(»BS(DtEL).LT.I.E-2)  go  to  700 

continue  - - 

nn,nC*C(l) 

IS  NN  active  at  this  time 


HVORu.aT 

HVDRO.rS 

hydro, v9 
hydro, SO 
hydro. St 
hydro, 5^ 
hydro, ss 
hydro. 5« 

HYDRO, S5 
-hydro. 96 
hydro, ST 
hydro, SP 
hydro, 59 

HYDRO. fcO 

hydro, 6l 
- hydro, 6? 
hydro. 65 
hydro, 69 
hydro. 6? 
hydro. 66 
hydro. nT 


IF fKlNDF(NN) .LT.R)  go  TO  lOS 

IF ( (T Im-tCHARCNn) ) ,tT,  l.t-i)  CO  TO  10b  

GO  YO  6«5 

continue  — _ 

mIsHIA(Nn) 

IE(mI.EQ.O)  CO  TO  110 

restore  aLU  variables  fur  event  NN  PREVIOUSLY  SAVEU 
THfSF  variables  aRe  RFCENERaTED  FOR  EACH  CaLL  TO  hYqrO 
ANd  need  not  be  saved  in  The  main  OVEhlAY,  

PSisPSIAfNN) 

VNrRHavNQBMAtNN) 

TZcTZ4(Nn)  — 

SRt8SRlA(NN) 

SI«SIA(Nn)  _ 

OIsT*OiSt*(Nn) 

H«BA(NN)  . 

mUmi  a(Nn) 

GO  TO  200  

continue 

get  distances  from  event  and  line  of  centers 


hydro, 68 
hydro, 69 

hydro. 70 
HYDRO, 71 

hydro, 72 

hydro, 75 
- HYDRO, 7tt 
hydro, 75 
hydro. 76 
HYDRO, 77 
hydro, 78 
hydro. 79 
-HYDRO, 80 
hydro, Bl 
HYDRO, 82 
hydro, b3 
hydro, 89 
hydro, 85 
-.hydro, 86 

hydro, 87 
HYDRO, 88 
hydro, 89 
hydro, 90 
HYO«C,9t 


HYDR0I9J 

CF| CULaTe  THf  OISTancFsDIST  FRO"  THF  LINE  OF  CENTERS  TO  THE  PDINT  mvdRO.bu 
ANq  the  OlStANCE»R*0  from  The  FINAL  EVENT  PUSITICN  TO  THE  POINT,  HYDRO, 95 


hydro, 96 

C»LL  RI8c*P(YIN(1,NN),xFv(1,nN),R)IYZ,0,SI»UIST,8b1)  hydro, 9T 

call  SuevtC(X£V(l»NNJ»PXYZiC) hydro. 98 
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210  2nn 

210 


227 
2J5 
255 
201 
255 
260 
26 « 

262 

265 
260 

266 
267 
27) 

275 

277 

506  650 


R60»X’'AGtC) 

c*lcui.*tf  Radius  of  evist  at  ti«e  cf  PASSACfc 

17  (8ltLT.-l,5j*REVfK))  CO  TU  700 

SRpaSRl-s) 

RR(fi&]N(MS)*SR0AR7v(NM)*Sl}/SRl  

R = A“AXlf(jEI^'t6'’^)»*PI7'l(P£''C7'N),R)) 

RHs01ST/r — 

RR«RA0/R 

PSj  IS  CaLCUlaTEO  tor  unit  RADIUS  AND  VELOCITY 
PSioO,5*rH*Rm*(1,-1./RR«»5) 

iPfPSI.Ci.l.RO)  GO  TO  050 

fIno  The  relative  time  

SRn»"SRD  

C*lL  SvZyr.y(9,VRZ(NN)»PSl.0iST,SR0»TZ) 

VNoRM«VRz(NN)/H 

nOpHALIZe  the  TIME  difference  TU  AN  tVENT  «ITM  unIT  RISE 

CONTINUE  

TRfLX=TZ.DTEl*VNORm 

aF,o<0  

YFbOIST/r 

ZFsSl  

C6l L ClPHF«(TRtL*.PSl,XO,ZO,NFL*C) 

iF(nFlaC.ED.I)  CO  TO  650  


00  Galilean  transformation. 


ZP.SRI-VrZ(NN)*DTELAZQ*R  — 

vF.XUAR 

C*|L  SuBvFC(P*YZ,XlN(  I ,NnJ  ,0l 

lF(XMAG(D)tLT.1.0)  0(11*100, 

C*l L LOCLAX(cvEC(l,NN),0,l,2,UVtC).  

NRA>i 

nCaRBs5  - ...  

NCt)ll 

POFdlaZF 

P0f(2)*Vf 

PCfCSIsXF  

*■1 

C*lL  mATmULTCIJVF  C.POF»P.nRA,nCARB,nCB) 

C‘lL  VEClINCa.XIsCI ,NNJ ,b,0,XOUJ 

continue 

iFfNUMX.EO.l)  CO  TO  700  

mUI 

s*ve  The  flom  field  parameters  for  event  nn 


mvdRC.rR 
HvnRi  , 1 01 

HVORU.lOl 
HVORO, lOi 

mvdho.iO] 

MVCPO, lOt 
HYDRO. lO' 
HviiBo.lOe 
HvrRO.tOl 
HVnRO.TOe 
hydro,  I0<i 
hydro. lie 

HYORU.m 

-hydro, 11 
hydro, 11 
hydro, 1 1 

HYDRO. 1 1 
HYcRu.iie 

hydro, 1 n 
hydRC,, 
hydro, ii<; 
hydro. 1 2( 

HYDRO, i21 

hydro. l2i 

HYDRO. l2i 
HYDRO, 1 
hydro, 1 21 

HYDRO. 1 ?t 
HYpto, i 21 
HYpRC.lZt 
hyuRI),  1 2<i 
hydro, 1 51 

hydro. 1 51 
►•Yo“0, 1 5« 
HYDRO. ill 

hydro. 1 5* 
hydro, 1 51 
hydro,  1 5e 
hydro, 1 51 
hydro, 1 5< 
HYDRO. l5) 

hydro. 1<I( 
hydro, 1«1 

-HYDRO. 1 Ui 
hydro, IRJ 
hydro,  I <l< 

HYfRL.lRl 

hydro, !«( 
hydro, 1 « I 
MYJRC.1«< 
hydro, 1 u< 

HYDRO, is; 
hydro, l5! 
hydro, i5i 
hydro, i5' 
hydro, l5< 
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HyDRfl 
Jl? 


ill 

SI* 

3t’ 

J?! 

iZ? 

J?b 

312 


6nn 




C-- 

339 

300 

303 

68S 

303 

5S1 

686 

SSI 

3S3 

C 

C 

C 

C 

3S5 

362 

67n 

362 

36T 

706 

36T 

C 



C 

c 

370 

373 

173 

37S 

377 

377 

001 

<10<1 

007 

«|3 

923 

906 

«25 

«2T 

030 

1006 

c 

c 

c . 

c 

c 

• 33 

•33 

a}0 

9 3*. 

Ml4(NN)lMt 

pSfl(MM)aPSI 

v^rs«»(^lg)»v^f|RM 

OIsTA(^^)«OIST 

rA(»,m)sr 

St»(N^)■7r 

COniTINuE 

TiACN^jPTOtlX 


SRlAtNMp^Rl-vOZCMNjADTEL 
jE  fA*I,Eo.l)  GO  TO  670 


S*VE  The  VECTnS  OlSPL*CE>'tNT  C*LCUL*TtO  fCR  EVENT  NNaSROU 


C*1L  XMlTfS.xrU.SPOUTl.D) 

CQ  TU  686 

COnTikuE 

C*lL  XMITC3,0,0,SPOUCI.I.)) 

continue 


iE(l«ct,o)  cq  to  too 


hycBC 

- HYpKC 

hyobc 

KCV?0 

mydno 

mYrho 

hydro 

--  HYO»C 

hydro 

HYDRO 

hydro 

hydro 

hydro 

- myoRC' 
hydro 
hydro 
hydro 
hydro 
hydro 

- HYDRO 

hydro 

HYDRO 


oc  the  Vector  addition  op  the  oisplacehents  for  pxyz  froh  HOLTiPitHYrRo 


events. 


C6(L  EvEnAD(SROU,NAI,PXYZ,XOU). 
COnTIHijE 

C*LL  XmIt(3,xou,XSaV) 

COnTihoE 

TImaTOBIo 


IS  TIH  An  event  TIME- 


iPfRPP.EO.l)  CO  TO  TOO 
TIhaTOBio 
nCdUsnCCuaI 
THpL  tHCru)«T0Rl0 
NCrT*sCCi)  ■ 

aUpxAc 1 n .a5*rmZe  RO 


C*|L  5CHCAtXru»».*LPH*»XS*¥) 
C»lL  X-TYt3,xSAv,XOU) 

C*Ll.  »NlT(i»XCIJ.XAR(lfX)) 
C*iL  *H|Tf 3.xrU.X*R»Y(I*X)) 
CUnTIhL’E 


jP  (NtAG(,^Aj).FB.K)N*I«N*I-l 

L»N*I  

continue 


T"f  node  S'lTfH  IS  set  TU  1 TO  INDICATE  THAT  The  PUInT  HAS 
undergone  St’HE  “OVEheNT,  * CHtCN_J3  .nadE._hERE  .TC -CETtHNjNE_ 

lE  tmE  hotiOh  is  Significant 


XN| 10,0 
“CoE«0 

DO  1200  |t«t.NCOT 

C*LL  SuPvCC(XARAY(t,n),P*TZ,Ci_ 


hydro 
hydro 
— hydro 
hydro 
hydro 
HYORU 
HYDRO 
hYD“C 

hydro 

hydro 

hydro 

hydro 

hydro 

hydro 

hydro 

hydro 

HYDRO 

hydro 
hydro 
hydro 
hydro 

HYDRO 

HYt'RU 

hydro 

hydro 

hydro 

HYDRO 

HYpRO 

HYDRO 

hydro 

hydro 

hydro 

I5YCBC 


ST 

SB 

S9 


61 

62 

63 

6H 

6b 

66 

67 

68 

69 

70 

71 

72 

73 
7a 
7b 
76 
71 

78 

79 

80 
FI 
82 
83 
80 
8b 
86 

87 

88 

89 

90 
Rl 

92 

93 
90 
9b 
9e 

97 

98 

99 

.200 
.201 
.202 
.203 
.200 
.20b 
• 206 
.207 
.208 
.209 
.210 
.21  I 
.212 


t I 


l?ftn  CRmTIMiE 


I''f*”J,CT.TR*vL)  MOOE»l 

PfcTUMN 

ENo 


MVCRO,?! J 

« 

HvnBo.ali 

hvdRU.JIT 


SUePPCCRAM  ieSCTH 


fusction  assignments 
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000000000000000000000000 


COITX 


SUrPQUTInE  EOITX(PxVZ) 

THIS  SURpbOCRAh  lOCPS  THRnUGh  THt  EVENTS  iN  STURACE  UP 
TO  THE  3THUL»TinN  Tint, 

THE  APRAv  OP  pertinent  events  iS  STCREO  IN  NCAQ 


INPUTS 


THprUGH  COMM  EVXOAT- 


Ef'lTX.z 
EtllTX,? 
ECITX.u 
EDITX.S 
EOITX.a 
EOlTx.r 
tOlTX.s 

EOITx.9 

EOlTX.lO 
EOlTX.il 
EOITx.i? 
EOITX.iS 
EOlTX.lfl 

EOITX.iR 

EOlTx.lfe 

- - f OlTX.lT 

TMJS  routine  IS  A FIRST  CUT  AT  EDITING  THE  EVENTS,  EOlTx.iS 

the  criteria  iiSTO  are  SIMPLE  AND  CAN  nEJECI  EVENTS  THAT  SHOULD  NCTEOlTx.iR 


NUMXaTOT*L  number  OF  EVENTS  UP  TO  TSIM 
tsih»simulaticn  time 

OUTPUTS 

NC/1GSARR4V  Oe  NaI  pertinent  event  indexes 
nAibnumBfr  of  Important  events 


BE  rejected  and  ACCEPT  EVENTS  THAT  SHOULD  BE  REJECTED. 


IF  an  event  is  ACCEPTED  AS  INFLUENCING  TmE  POINT  HlSTORT»IT 
SHOULD  RE  checked  after  THE  OETAlLLD  CALCULATIONS  HAVE  BEGUN  IN 
subroutine  HycRO. 


EDITx,?0 

EOITX.pl 

EDITX.?? 

EOITx.pJ 

EDlTX.pU 

EDiTx.pS 

CEClTo.p 

CEOTO.s 


time  dependent  model  parameters 

COMMON  /CEOIn/  NF.  InOXFISU).  RTFCSU)*  RLEC50).  HFt50)»  uCFISO). 

1 GlF(SO).  h'*XF{SO)»  hminF(SO).  KlNCFlSoJr  T I LTF  (bO  J , r.f.cjTG  . « 

2 AGEISO),  NBTA,  INDXO(IOO),  DLAbL(lOO),  hCRcIOO),  uEQTO.S 

5 MDR(IOP).  RTbSdOOJf  RLBS(IOO),  hRS(100)»  HN  fc  S ( 1 0 0 ) , GEH  T J , 

R CCETA(lOOJ,GLBTAClOO),IF(,;0),TCMAR(50),MR610(5g),  GEOTO.T 

5 XFRC5n).»FSI50).ZFR(S0).RUTlS0) 

Cai.MON/EvXO*T/NUHX*NAl,X*R(J,  I0)f  nCAGCIO)  » T5im,nCOT*TmBU(10) 


C0MHON/RiD0AT/XIN(i,10),XfV(j,10),RElN(l0),RtV(l8),VXTI{3»10), 
1 VRZtiO)»CVECfJflO) 

COmhON/GeOHD/Tme TA.RR.RH, VN0RM,0TEL,XR0X,TZ 


EVENTX  P»PAMETERS 

COmmUF  /eveNtx/  NX,  IDX,  TR(SO),  hbCSO),  GCfltSO),  GLB(50), 

1 iDCADfSO),  RhOHCSO).  mSR(SO),  TEmB(So),  V R I St f bO ) , I VE NT  X . u 


GEOTD.A 
tVxPAT.Z 
tVxOAT . 1 
hadDaT.Z 
raddat, j 
NApDAT.U 
CEOHO.p 
GEOMO.  J 
EVENTX. 2 
EVENTX.  J 


-2 

S 

a 

» 

T 

a2- 


2 ROZEROIbO),  PHZEMO(bO),  RUZEROISO),  BXB(bO), 

S BTH(SO),  fcZB(bO)E  LMVbCbO),XALPHA(bO),  hALCH 

COmhON/RfpQRT/NAFECT  

dImENSITn  px¥7(3),C(J),R(S),U{J),HIJ) 

N»HtLTST/F0lT3/PXTZ.C,L,SRl2  . . 

begin  EVinT  LOOP 

NCn»0 

nAeFCTeO 

NAf  »1 

NCtG(|)R0 

DO  1000  l*I«NX 


SHl2«*MAc(CVFCtl,L)) 

- ir(SRl2.LT.l.ERJ_C0  _T0.1000_ 


EVENTX. R 
EVENTX. 6 
_NOv20A.71 
EOITX.JI 
EDITX.jZ 
EOITx.tJ 
tPiTl. JR 
ECTTx.iS 
.NC'V?!)*.72 
EDI Tx , J* 
eoitx.tt 
EOlTx, J8 
tPITX. J9 
EDlTx.dO 
.EDiTx.ul 
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--  18  THIS  fVfcHT  *CT1VC  I,  E,  NOT  * HERtifcO- EVENT - 

ft  Continue  

LU*D  event  CiRAMETEPs  F0i»  the  EVENT  L 


nIil  This  l-event  ever  m*ve-*n  effect-on  point-civenj 

begin  CiirULlTION 

C»LCIjL4.Te  distance  of  CLOSEST  APPROACH  

C*lL  0IScAP(XINfl,L),XEV{l,L),PXYZ(l),H,Sl,0IST,SR2I)  - 

RDtS  »OIsT 

T0fS»(SO?1-Sl T/VRZ(L) 

^lE  (SI  ,LT..l,S»fif  V(LnT0IS  = ee8fl, 
iFf  (SRZl.sn.l  T.-l,S«REV(t))T0ISse55S, 

THf  RADIUS  RfV(L)  should  NUT  BE  USED  HERE,  R The  RADIUS  CF  THE 
EVpsT  nEa»  IhF  point  of  ISTERFST  ShUULU  faf  USED.  HOWEVER#  THIS 
RArjUS  H REOUpBES  extensive  CAlCUlXIiON  NUl  uEEmEO  JUSIiMEOAT 
ThtS  point  In  The  mIsTCRY  PROCESS, 


HH,RDIS/REV(l)  

IF  fRH.GT.I .81  CO  TO  lOOO 

iFfTDiS.r.T,  7T77.)  Go  TO  1000  

TlpSTtTSlM.TBfLlAl.S.REVlD/VRZlL) 
IF (TOIs.CT.TTEST)  go  to  1000 
NCpiNCC* 1 

nCaC(NCC)«L  

NAjinCO 

continue 

NApFCTaNAt  

return 


EDlTx.aP 
EDITt.oS 
EDITx.uR 
tniTx.uS 
EPlTx.ub 
EDITx.oT 
FDITx.ue 
-EDiTx.aR 
EDI Tx ,50 
EOITX.5I 
EOITX.S? 
EOITx.sJ 
EDITX.S# 

- EOITX.55 
NOVIS.I 
EDITX.sT 
EPlTX.fcl 
EOlTx.fc? 
fOITX.bJ 
EOlTX.Ntt 
EOITx.fcS 
EDITX,66 
EOITx.67 
LDlTx.fc* 
EOrx.sR 

EOITX.tO 

E0lTx,7l 

EOITX,t2 

E0ITX,7S 

EOITX.t# 

ED1TX,75 

E0iTX,76 

EDITX.7T 

NOv20A,7J 

EOITx.xe 

EDITx.tO 
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i 

1, 

J 


oori  or>r»  nnooo  onno  rir%on  of%r»  oo 


s>jbpuut!me  EoiT*(PxyzJ 

This  S'jHpPOr.RAH  LOCPi  THHroGH  THt  tVtHTS  IH  STCBAGt  UP 
TO  TriF  SlMULATins  Tint, 

THf  ABRAy  OP  PtRTlStsT  EVENTS  iS  STOHfU  IN  NCAC 

inputs 

TMrpuGh  CDmh  EVXOAT 

nUhX=toT4l  number  rr  events  up  to  tsih 

TSiMsSiHui ATlON  TIhE 

outputs 

NCiGSAMRiV  Ce  NaI  pertinent  tW£NT  iNOtXfS 
nAiSNUmBer  CE  important  events 


ET'lTx,? 

EIITX.X 

EfTTA.a 

EniTx.s 

EGITx.n 

EITITa.t 

EOiTx.p 

EOlTx.p 

EPlTx.iO 

tniTx.it 

LOITx.i? 

tniTx.i 5 

ET’pX.iU 

tniTx.is 

t r I T X , t N 
EUITx.i T 
tOTTx.iB 


T'MTXtl  T 

TMiS  RCUTTNE  IS  A FIRST  CUT  »T  FUITINU  ThE  EVEMS,  tOITx.i? 

TNf  criteria  iiSFR  are  SiHPEt  AnO  CAN  HEJtCl  EVENTS  THAT  SHOOED  NfTEDlTx.ig 
BE  rejected  anU  accept  events  That  ShOuuO  at  REJECTED,  EDITx.pr 

EDITX.Pt 

IE  AN  event  is  accepted  as  INELUEnCInG  ThE  PUInT  MISTORV,!!  fPTTl.p? 

SHnOLO  BE  CHECXeO  AETFR  the  OFTaILEU  CALCULAI IONS  HAVE  BECoN  IN  ECHx.?? 
SUijROUTInE  HypRn,  EP'ITx.pU 

EDI’x.pS 
SE' Tu,2 


TIhE  dependent  hhDEL  PARAMETERS  GE'Tu.j 

COmmUN  /r.EOTn/  me,  InOXEiSu),  KTE(SU),  PLEiSJ},  hE(SO)#  mCEISO).  CEtTo.x 

1 GlTISO),  m AxE(SO)#  HHitFCjO),  XTNrMSO),  I I L T v ( SO  ) , l-E  ( T „ 

2 ANttSD),  ntTa,  1NL>XD(10C),  ULA'.iLhoOj,  hl'RiIOC),  mEi'Tj.n 

3 H.TRdflD),  RTBSdOf'),  ROasdOd),  MPxdCO),  RNr  S (1  00  ) # r.F  T.,  , M 

“ GCHT4ttnO),GLE'TA(ll'0),IA(sJ),TCi-*H(^U),MNOlClSO),  Ofr.Td.T 

C tfAUrCA^.VCiarC.t3.7L^rC,.fl3.M..TEt.ET3  D-r  .Vr.  . 


S XER{5n).vfRlR0),IER(Sy),HuT(S0) 

COmmON/E  vxO»T/NUvX,nai,  X*R(  J,  10)  f nOA  Jdo) , ISim.nCDT,  THHUdO) 

COMHaN/R»OOAT,/XIN(  j,l0),XE  VCi,10),REIN(10},BtV(t0),VXYI(i,lO), 
1 VRZdn),CVEC(3,tO) 


H ADL'AT  a 

COMH'JN/GErMD/TMETA.RR.HH,  VNPRh.oTEL,  «PuX,TI  GE'.mu.P 

GEO'-O.  ( 

eventx  Parameters  EvANtx,2 

COmMUT  /event*/  NX,  lOX,  THIS"),  mbCSO),  GCBIS"),  G0b(5"),  tvEMxIs 

1 iDGADISO),  HmHH(SO).  mSP(SO),  TEmPISu),  vRISltSO.LvE* TX.u 

2 RCEEROISV),  Rh/eho(S0),  HU/tROISO),  HXcjTSO),  tvE‘‘T«.S 

3 PTHTSO),  PEHISO)*  LhvI'ISU),X«oFmA(SO),  KALCH  tVENTX.*. 


GF'’T';,p 

EVeDaT 

tVxi'AT.  t 

hadT'AT,? 
baddat  . j 

HADl'AT.a 
Gtf.MU,? 
GEO'-O.  ( 

EVANt*_2 

tvtMx.3 


COM“riN^HEnnBT,EAEFCT 

dImENSITn  PXy7(d,c(  j)  ,«(  3),r.(3)  ,HtS) 
N»MtLrST/EDiTX/PxVI,CiL»SRl? 

begin  event  loop 

NCn»o 
nAeE CTsO 
* A j*l 

NC*cd  )sn 
po  1000  I *1  ,*JX 

5H,/*X>»r.fcvTrd,L)) 

lP(SRta,LT,t,t«)  GO  TO  10"0 


tVENTX.*. 
nDv/u A , 7 1 
tDlTl,d 
tOTTi, ,P 
f DITx, d 
LrjTi, j« 
EDTTx , tN 
^nv?;A.7^ 
tOITx.x* 
IDITx,  tT 
t"! Tt. d 
EDlTl, to 
IDITx.un 
t"ItX,al 


L 


Syzvcv 


SUbBOUTInE  5Y?VGYfH,VB,PSl,OIST,SBO,TOUT) 

THIS  SU?PBOCil«M  c'LC'iLATtS  THf  0«IFT  FUNCTION  a TOUT 
FOR  * SIvEN  POINT  relative  TO  A SPHERE, 

inputs 

RBrAOIUS  of  current  event  AT  TIHE  OF  PASSAGE  OF  TP 

vRsRISE  velocity  - 

PSibSTreaM  functio  V*LUE 
OUTPUTS 

TZaTIHE  cCNSTaNT  CONTOUR  RELATIVE  TO  EVENT  CENTER  OF  ZERO 
COmpON/CONST/RE,PI,ThRO,PZEPO, AZERC,UZERO,T2ERO,FI02 
COmhON/CeOHO/THETAjRRjHH, VNORP,nTEU,XROX, TZ 


0*T*  PI02/I.S707B6J/ 


OlsT  IS  A POSITIVE  OEFINITE  QUANTITY  BY-.Tht.  OEFIKITICN.  C 

OP  THE  COORDINATE  FHaNE 

pStbAMAY) fPST.s.f.T) 

oIstaahaxkoist.i,)  -- 

TM7TAsaTaN2(OISTiS'RD) 

IFlAGbO  

VNnRHBVB/R 

ROrSURT(2,0»PSI)  --  — 

COsECs1./SIN(ThET*) 

COtANbsRo/OIsT  - — 

BOBaHO 

RHaOlST  

RRsCI3T*cOSEc/R 

PRINT  IObO,Bb,R,ROA,sRO.PSI 

FORPATdK.aWRs  AlPEia.b,*  hs»IpEl2,5»»  RU**SRD,PSia  *1P3E12,S) 
IP(RO.LT.I.S) 

l*MrZ®tU<itt<luu*f  5.*Rn«PP)«  ALOGTu  .5bR/Ra)-2,60UY0,nR0T»RO*RO 

iPtRO.GE.l  ,5)XROZB.<‘417''/RC*»»5“,'»IR/RU***OYI,9  026/RO*,»»H 

lFfTHF.TA_CT,2>17, OB. Theta, LT.O.bZibB)  CO  TO  200 
iPfBRiCTll .75)  GO  TO  50 

uTzeAB0Z/2,*.Nbho6ft*(l ,*CBUA»»2)/1,)»AL0C{ tAN( theta/2,) ) 

J •,22222»R(!A»»2»C0SIC*CoTAn  

1 ♦.2222»B0A»»2»CnSEC»C0lAN 

TZaUTZ/VNORM 

TQt«TZ 

iPl*'**!  - — — - 

IF(BW.GT.I,5)  GO  TO  tOO 

iFLAGeo  . , 

GQ  TO  30o 

continue. 

continue 

UT7a-RCA«C0TANTXR0Z  - ... 

TZbUTZ/VnORh 

iPflPLAG^FO.n  CO  TO  500  _ 

GC  to  300 

continue 


SYZYCY.2 

SY/VCY,! 

SY7YGY.<| 

SY2YGY.5 

SV7YGY.fr 

SYZYGY.7 

SYZYGY.S 

SYZYGY.R 

SYZYGY.tO 

SYZYCY.l 1 

8YZY6Y. J2 

SYZYGY.li 

CONST. 2 

CONST.! 

CEpHO.2 

GEpHO,  j 

SYZYGY.lfr 

SYZYGY.IT 

SYZYGY.I6 

SYZYGY.IR 

SYZYGY.20 

SYZYCY.21 

5YZYGY.22 

SYZYGY,2i 

SYZYCY.2R 

SYZYGY.Z'i 

SYZYCY.26 

SYZYGY,27 

3YZYCY.28 

SYZYGY.29 

SYZYGY.SO 

SYZYGY.Tl 

SYZYl,Y,S2 

SYZYGY.53 

SYZYCY.3R 

•'vr.Y.35 

'.37 

SY»  GY.3B 

SYZY6Y.S9 

NOV20.2 

NOv20.! 

SYZYGY.ai 

SYZYGY.H2 

S»ZYGY.U3 

SYZYCY.iiR 

SYZYGY.ab 

SYZYGY.Rfr 

SYZYCY.gZ 

SYZYGY.ua 

SYZYGY.UR 

SYZY6Y.50 

SYZYGY.Sl 

8YZYGY.S2 

SYZYGY.53 

.SYZYGY.S# 
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syzygy 

?n(, 


?53 

?5«- 

25S 

261 

266 

266 

270 
271 

271 
271 
277 

300 

301 

302- 


RRc.RR 

- B»p»ABS(pP-H  — 

UT7i-BR*.i^(,^;j,»(PO*RP),,2/(«R*»i.l  ,3.,1  1 t Ul»(3,*WO» 

1 *LCG(P48/SCHT(PP«BBaRBa1)  )•,  192<iS« 

2 /(1.a2.*HP)} 

pp,.pp  

uT7»UTZ»r 

YZsUt7/Vp 

I*"  fTMl  IA.lT.piCZ)  Tz*-rz 

I7  fTMtT*.LT.Pl02)TZsXR0Z/VN0R*«ATZ 

loo  continue 

TZsTZ»VN0RM  --  

TOuTBTZ 

RtTU86  

son  cOictimue 

TZ=((Rh.j  ,5)*TOT*(1.75-RR)*TZ)»<li 

TZsTZ»VNoRM 

TOuT=TZ  

rEtURN 

ENf> 


• 2)* 

1.7320% 


SYZVCY,S% 
S»Z»GY.S6 
S^ZYGY.SZ 
SvZYGY.Sd 
JvZ'tY.SR 
SVZ''C».60 
S»ZYGY,6l 
-SYZVt.V.sZ 
NOvU.i 
5YZYCV.6% 
60v2o , a 
8YZVr,Y.66 
SYZYr.r.o7 
- 9YZYCY.68 
SYZYCY.69 
SYZYGY.70 
SYZYGY,72 
SYZYGY,71 
SYZY6Y.74 


SU8RB06PAH  lEI'CTM-  . 
00012 


H 
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fi  no  noonoononnoon 


ClPHEp 


SUpROOTInE  ClPMFR(T,pP,X,y,KEL*C) 

TH,S  SiJPPO'iTjAt  C*LCiiL*TlS  ToF  PCSHIUN  C'f  * PUt».T  IN 

COrPOIf^lTFS  tiXfO  “ITh  The  SPMt»t  FOk  a given  STPE*** 
function  PP  and  a given  U»IFT  FUNCTION  T ISEC) 

inputs 

TpoSI^T  function  In  SEC  GCAIEO  TO  A UNIT  «tSE  yFLOClTV,— 
pasTREA**  function  FOR  UNIT  EVENT  RADIUS  ANy  UNIT  VELOCITI 


CIPHER 

CIP-E* 

ClP-E* 

CIPHER 

cipher 

ClP"tR 

CIPHER 

cipher 

CIP-ER 

CIP-ER 

CIPHER 

cipher 

CIPHER 

cipher 

CIPHER 


.f 

.T 

.S 

.T 

.10 

.11 

.!« 

.!!> 

.IP 


OUTPUTS 

*Sj.»3QHpi^,ATE  relative  to  EVEhT  CENTER  AND  PERPENDICULAR  TC 
RISE  VELrrlTV 

r=LCNCiTuf)iNAL  coordinate  at  time  t in  uniis  Relative- to — 

UNIT  RADp^S  AND  RISE  VEIOCITY. 

NELlOa^LAC  Tfi  INDICATE  -hCIhER  anT  Cm*nGE  in  POSITION  H*S  CCC'.iKREOC  I phe  R . 1 7 

CIPHER  18 

COhpON/'ConST,re,PI,TmRD,PEERO,AZERO,DZERU»IZERO,PI02  - CONST, 2 

const, T 

COmMON/GeOhO/THF  TA,RR,RM,  VNORM,OTEU,XHOX,  TZ GEPH-O.Z 

CFO-O, 


c 

. 



CIPHER. 21 

dimension  Titjsi.psds), 

CIPHER. 22 

0*T*  TI/.IO, 

l«6f-UR9*l 

t 9*  i**# 

cipher. 23 

A 

, ___  M 

uO,,bO#/  - 

rTOHFfl  PN 

0*TA  PS/S.l- 

cr##C 

7. 1 .?5E-S,b. 

E-b, 

•b«t»S9l9l2SF»29 

CIPHER. 26 

A 

? ttm?0  • 

04ib«  •0H#09  \ 9 160# 

,.5,1,125,2 

.01/ 

CIPHER. 27 

data  NP/tS/, 

NT/58/ 

CIPHEP.28 

OAyA  TlAeT/h 

1 8. /, PL*S3/1 

. 01002/, RlZ/IOO, 0660/ 

— 

CIPHER. 20 

OAtA(Xx(ij,1 

*1,11“)/ 

CIPHER. 30 

R 

,D,00?T8» 

5.01770,  . 

-0,02677, 

3.05222,— 

2.57502,  - 

-2.11066, cipher. 31 

R 

1,03880# 

1.61102, 

1,30001, 

1.23051, 

1.0/335, 

1,05627, CIPHER. 12 

. - R 

l.03?(i?. 

1.01830, 

1,01021, 

1.00312, 

1.00172, 

1.00070. cipher. 33 

t 

1 ,00028, 

1.00000, 

1.00002, 

1.00001, 

1.00001, 

1,00001  , CIPHER. 30 

p 

1 .0000 1 # 

1.00537, 

1 ,50866, 

5,30736, 

7.38‘»30, 

0,ifl586,ClPHER.35 

R 

13.38207, 

16.39200, 

21,38127, 

31.38060. 

01 , 180O7, 

10.00278, CIPHER. 16 

D 

5.01770, 

0.02878, 

3.05223,  - 

2, 57502, 

2.11067, 

1.03890, CIPHER, 37 

R 

1.61103, 

1.30001, 

1.23051, 

1.07330, 

1,05626, 

1,01201, CIPHER. 38 

R 

1 ,0l»TO, 

1.01021, 

1.00312, 

1.00172, 

1,00070, 

1.00C2O, cipher. 30 

R 

1 .OOnflO, 

1.00002, 

1,00001, 

1.00002, 

1 .00002, 

l.OOeiu, cipher. 00 

R 

1.70788, 

3.56008, 

7.52052, 

0.52620, 

11.52005, 

15.52276,LIPHIR,01 

R 

i8.522iU, 

23.52150, 

33,52113, 

03.52U05, 

10.00278, 

5,01 771, cipher, 02 

R 

R, 02870, 

3.05220, 

2.57500,  . 

2,11068, 

1.03882, 

1.61 105, CIPHER. 03 

Q 

1 . i<iOo3, 

1.23053, 

1. 07337, 

1.05620, 

1,03203, 

1. 01832, CTphlR. 00 

R 

1.01023, 

1.00310, 

1 .00173, 

1 ,00072, 

1.00030, 

1,00010, cipher. 05 

Q 

l.OPOfl'i, 

l.OOOOo, 

1.00030, 

1 ,00650, 

1 . 1 0602, 

2.52o0u,cIphER,06 

R 

II. 06082, 

e.o5lU8,. 

10.01.006, 

12,ou77o, 

16.00632, 

10, 00580, CIPHER, 07 

R 

20«<Hi5ti  , 

30.00080,' 

00,00073, 

1 0, cn20o. 

5.01/05, 

O.02Or8.ClPHER,i;8 

R 

3,05262, 

2. 51500, 

2,11522,  . 

1.03030, 

1.61260,  - 

J,30073/C1P-EO.OO 

0*T*(«xf Jl.l 

SI  15,228)/ 

CIPHER, 50 

R 

1.23522. 

I.OT5O5, 

1.05680, 

1.032OO, 

1.01870, 

1.01068. cipher. 51 

R 

1,003S7, 

1.00217, 

1.00118, 

1 ,00066, 

1,00101, 

1 ,00268,CIPHE  R.52 

p 

1.01030, 

l.l“02n. 

1 ,76002, 

2.65773, 

0.6C878, 

6, 50656. CIPHER. 53 

R 

1 0 »580uS, 

12.58812, 

10,58730, 

18,58638, 

21,58600, 

26,58562, CIPHER. 5o 

R 

]6t5fl528. 

06.58513, 

10,00328, 

5,UI 860, — 

0,03008, -- 

- 3, 05383. CIPHER, 55 

50 


cipher 


R 

2.116O0, 

1 .001 lO, 

1 .61060, 

1 . 30?00, 

1, 23738. ClP*-tW. 5b 

R 

1.07578, 

1.05650, 

1.03050, 

l.u?u3a, 

1,01213. 

1 ,00507.ClRi“tR.57 

P 

1.0038?, 

1.00335, 

1.00032. 

1.00009, 

1.03037, 

l.l?317,ClP-tK.56 

R 

1.80081 , 

2.57511, 

3,53763, 

5.51328. 

7.S05O5, 

ll, 50000, CiPHtR. 50 

R 

iS.RO^RO. 

15.00828, 

10,00705, 

?2. 00700, 

27,00672, 

37,UR«.37,ClP6tP.60 

Q 

«7.U()6;>n, 

1 0. 003On, 

5.01003, 

o.0315o. 

3.0558O, 

2,57025,ClPHte.9i 

R 

P.llOfcO, 

1 .00018, 

1.61803, 

1 ■ 30651. 

1 .20p6e. 

1 ,0788o,ClP6EH.62 

— - R 

1.0615?, 

1.03727, 

1.02201, 

1.01072, 

1.00811, 

1 .00750, CIPHER. 63 

R 

i.oooio. 

1.01535, 

1.03023, 

1.13310, 

1 . 35890, 

2.1  3803.C!P6f>;.6O 

R 

3l  07bht  9 

U , 0*5  1 7<i » 

6.03305, 

8.02707, 

1?, 02270. 

lo ,02165, ciPHtP. 65 

R 

|b.0?0O5, 

?0.0?01  1, 

23.01073, 

26. 01 0 lo. 

3=. 01*00, 

«6, 018/5. CiPHtP. 66 

R 

10.00U78, 

5.02167, 

0, 03370, 

i t oSbbb  # 

2.58255, 

2. 12550,ClP6tP.67 

R 

1 .RU8-\7, 

1.62277, 

1.35160, 

1 .20000. 

1.08318, 

1 .06571  ,cIphK)_60 

p 

l.nai?! , 

1 • 02b7u# 

1.01866, 

- 1.01333, 

1.01030, 

I,02l3o/CIPH18.60 

D*T*C)ixtt),  I 

8220,3021/ 

CTPHtB.70 

_ - - R 

l.OSRuO, 

1 ,10053, 

1 ,26610, 

1,60857, 

2,07107, 

3,0?765,ClPHL«.7l 

R 

<I,<10‘»37, 

6.3O035, 

8.36807, 

12.38385. 

1 0. 39?7e, 

lb,382C5.CtPNto.72 

...  R 

20,38115, 

23.38070, 

28,36020, 

38,37083, 

08.37050, 

10. 00727. CIPHER, 73 

R 

5.02662, 

0,03080, 

i # 06feb  ^ # 

2.50103, 

2.13072, 

1 .06028,C Iph£R.7o 

R 

!*636?7# 

1.36608, 

1 .26002, 

. 1.00502. 

1,0/808. 

1,05311, cipher. 75 

R 

1.03877, 

1 .03l7i , 

1.03332, 

1.00205. 

1 .07260, 

1.13oO<»,cipHEB.76 

R 

1.20180, 

1.60000, 

2.01560, 

2.03661, 

3.00663, 

0,80223, CIPHER. 77 

R 

6 • fl7R  j 6 # 

6,87506, 

12.86=55, 

10.66732. 

16,80605. 

20. 86533, CIPHER. 78 

R 

83,86070, 

26,86010, 

38,86352, 

00.8631b. 

10.01077, 

5. 03355. CIPHER, 70 

R 

«,0«Pu3, 

3.07786, 

2.60501, 

2.15021. 

1.07687, 

1 ,65507. CIPHER. 80 

„ R 

1.38633, 

1.28071, 

t. 1 louo, 

- 1.00627. 

1,07 1 10, 

1, 05772. CIPHEM. 81 

R 

1,05320, 

1.06016, 

1,002??, 

1.15817, 

1 .27006, 

1 ,o675o,ClPHfP.*2 

R 

1.87065, 

2.30020, 

3.20878, 

0.22308. 

5.20070, 

7, 10660. CIPHER. 83 

R 

0,IOOu3, 

1 3, l8a7o, 

15.18323, 

17,18210. 

21.18068, 

2o , 1 7006 .CIPHER .80 

R 

?0,170,?, 

30,17616, 

00, 1 7763. 

1 0. U 1 526, 

5.0O2O5, 

O.OSOoo. CIPHER. 85 

R 

3.0O??0, 

2, 62170, 

2,16000, 

1 .00800, 

1 ,670o5, 

1.0122<>,ClPHEa.86 

R 

1.30688, 

1.13016, 

1.12080, 

- 1.00603. 

1 .08055, 

1 .0aoOb,CIpHER.e7 

R 

1 . t20?8. 

1.16120, 

1.261)03. 

1.01020, 

1,67608, 

2,0S02o,ClPHER.ea 

R 

2,5263?, 

3.07281, 

0.00806, 

5.03002. 

7. 02030, 

0,oi327/ClPHtR.S0 

0*T*C«X(t),l 

=505.056)/ 

CIphEH.ro 

R 

1 3. 0o6u0, 

15.00055, 

17.00316, 

21.00126. 

20,00020. 

20, 30016. CIPHER. 01 

R 

30.30783, 

oO. 30707. 

10^02075, 

5.05330, 

0,07287, 

3,10063, CIPHER. 02 

R 

- 2.60205, 

2.10306, 

2.02367, 

_ l»7060bi 

1.00385, 

1,33601, CIPHER. 03 

R 

1.17007, 

1.15211, 

1,12820, 

1,11055, 

1,12006. 

1.18233. CIPHER, Oo 

R 

1.23023, 

1.36510, 

1.53785, 

1.82078, 

2.20005, 

2, 60178, CIPHER. 05 

R 

3.60010, 

0.61068, 

5.60000, 

7,58o23, 

6,57601 , 

13. 56773, CIPHER. 06 

R 

15.5653T, 

17.56361, 

21,56116, 

20,5500o, 

20.55800. 

30. 55661, CIPHER. 07 

R 

00,55557, 

1 O.fl347o, 

5,0808?, 

0, lObOo, 

3.15560. 

2.60il3,ciPHtR.oe 

. - - ..  ......  R 

2«2SUA0f 

2.0»778, 

1.78056, 

1.52316, 

1.02013, 

1.25288, CIPHER. 00 

R 

1,235(0, 

1.2I020, 

1 .21?08, 

1.22862, 

1.32137, 

1,5ro35, CIPHER. 100 

R 

1,55053, 

1,70818, 

2,05077, 

2,07381, 

2.02550, 

3. 86007, CIPHER. 101 

R 

0,80030, 

5.82216, 

7,80101, 

O,7O0Po, 

13.77701 , 

15, 77o3o, cipher. 102 

R 

17,77163, 

21.76770, 

20,76577, 

20,76335, 

30.76030, 

00, 75866, CIPHER. 103 

R 

10.05?6?, 

5.1 I50R, 

0, 15013, 

3.20020, 

2.75753, 

2, i?0?3, CIPHER. 100 

, , > R 

2.16761, 

l.e70«6. 

1.62200, 

1.52265, 

1.36017, 

1 ,30306, cipher. 105 

R 

1,32585, 

1.32020, 

1.35553, 

1.06565. 

1.55250, 

1,71/1 1 .CIPHER. 1 06 

R 

1.'’1557, 

2.21870, 

2.6  181?, 

3.08065, 

0.02010, 

0,08  367, CIPHER. 107 

R 

5.060(10, 

7.03301, 

0,01738, 

13.00030. 

15.80520, 

17, 80126, CIPHER, 108 

R 

21.88563, 

20.88260, 

20,87005, 

30, 87050, 

00,87106, 

I0.1IO5//CIPHER.10O 

0*T‘l»Xt J)/I 

=057,525)/ 

ClPMEO.nO 

R 

5.23628. 

- O.2O710,- 

_ 3.30586, 

2, 07102,.. 

_ 2.57531,_ 

_-.2,<l28«5»ClPMtP.lll 
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C1P«F8 


2>l6lS?f 

1.94J JA, 

1.65394* 

1.7ll7i, 

1.69802, 

1 ,68b04'*ClPHtR 

.112 

R 

1, 69475* 

1.72397* 

1.69103,- 

1 i92b3o.- 

2.06517, - 

2, 26605, CIPHER 

.11* 

R 

2*55tS5. 

2.9449t,* 

3.36831* 

4,26769, 

5.20623, 

6, 16499, CIPHER 

.11“ 

R 

». 1)373* 

1 0.ft8311  * 

19,04899, 

16,03636, 

18,03018, 

22. 01812, CIPHER 

.115 

R 

^S.nnoB* 

30.n042n, 

39,99458, 

49,98864, 

16,20067, 

5, 40025, CIPHER 

.116 

R 

U,U9U9S* 

3.69163* 

3.29767, 

2.68767, 

2.75615, 

2. 52048, cipher 

.117 

R 

2.32974, 

2,?5454* 

2.13509, 

2.12381, 

2.1 1968. 

2. 12328, CIPHER 

. 1 18 

R 

- 2. Iu9a7*  - 

2.25226* 

2.  32699,— 

2,46976,- 

- 2.62943,  - 

2.b8317,ClPHEO 

.119 

N 

3.2U26S* 

3.63h23* 

4,48906, 

5,39916, 

6.32667* 

6,2450) .CIPHER 

.120 

R 

to. 194po, 

t«. 13592* 

16,11770, 

16.10355, 

22,06299, 

25, 07190, CIPHER 

.121 

R 

JO. 058117, 

90.09150* 

50.03190/ 

CIPHER 

.122 

— - 



CIPHER 

.'2* 

•7 

11 

11 

12 

12 

16 

21 

23 

26 

27 

31 

32 

33 


16 


20ft 

21fl 

C 

c 

c 


BEcIN  IKTERPoLaTION  or  X in  T*BLE  *T  time  I *No  stream  PSI, 
lNpLG»0 

pspp 

COmTUUE 
M0pF«0 

iifp.LT.psd))  Norr*! 

ir(KOfE,EO,J)PePS(l) 

PO  200  M32»mp  

m2,n 

irfP.Lt.PS(N))  CO  TO  210  

COmTInuE 

NUpFa2 

continue 

FP«(PS(N?)-P)/(PS(N2)-PS(N2-1)) 

TMpBE  need  Be  no  PESTBICTIONS  on 
OtiTION  IS  LOCABITmmiC, 


FP  IF  TMt  Interpolation  • 


CIPHER 

CIPHER 

CIPHER 
- - CIPHER 
CIPHER 

cipher 
cipher 
cipher 

CIPHER 

CIPHER 

CIPHER 

CIPHER 

CIPHER 

CIPHEN 

CIPHER 

cipher 

CIPHER 
EXTRAp.C IPHE  R 
CIPHER 


37 

ITqFFbO 

cipher 

uu 

37 

lF(T,LT.Tl(l))ITCrFB(  . . . .. 

cipher 

95 

93 

DC  300  ib2,Nt 

CIPHER 

ub 

95 

Mlgt  - . . -- 

CIPHER 

47 

ab 

Iff  (T-TKDJ.lE.I.E-o)  go  to  3l0 

CIPHER 

aa 

S2 

30ft 

rC.jTiNiiE 

ClRHtR 

49 

59 

lTnFFB2 

cipher 

50 

55 

3tft 

continue  . . 

cipher 

51 

55 

lPAO»NOFFtl 

CIPHER 

52 

57 

N*OOSITCFF»l  

cipher 

53 

61 

GC  to  (9tO,92n,990),N*00 

CIPHER 

59 

67 

91ft 

GU  TO  (960,960,500) , JPAD 

CIPHER 

55 

76 

U2ft 

GO  TO  (560,560,500), IPAO 

CIPHER 

56 

105 

99n 

GO  TO  (520,520*500) , iP*0  

cipher 

57 

119 

96n 

continue 

CIPHER 

58 

122 

FT8(T1(Ni)-T)/(TI(n1)-7I(NI-1))  . 

CIPHER 

59 

126 

FTaAHAxi(0,,AMINl(FT,l,)) 

cipher 

bO 

135. 

Al,FP*FT 

ClPhtW 

6l 

136 

A2x(l,"FT)»rp 

CIPHER 

62 

137 

A3BFt*(I..FP)  

CIPHER 

63 

191 

A9,(1,-Ft)«(,.-Fp) 

CIPHEP 

69 

192 

XlaXXfNl.l ,N2>1)  . 

CIPHER 

65 

199 

X2aXX(Nl,N2-]) 

CIPHER 

66 

195- 



x3aXX(Nl>l,N2) 

CIPHER 

67 

2R 

25 

26 

27 

28 
26 

30 

31 

32 

33 
3« 

35 

36 

37 

38 

39 

40 
4| 
1*2 
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1 


to7 

ISI 

»7« 

171 

172 
»7T 

2e« 

2fl* 

2tl 

220 

222 

225 

225 

2S0 

231 

232 
232 
230 
230 
237 
200 


200 

2St 

252 

252 

253 
250 
250 
255 
261 
261 
263 
270 
276 
300 
312 
312 


«7n 


oha 

C 

C 

c 


500 


520 


530 

C 

C 

c 


tl< 


«0»*2))- 


*0»«X(Ni 
BC»S'<OT(oPi(ii< 

R2sOU*Ha 

R5eOO*«2  - 

«TF»2s2..P»Hj/(R3»n 

lii^ORT  t*TE>'25  - 

vTirMiEwp.xttM? 

1> t.o.O)YTEh?«0,0 

V5.sr)RT(vTl>V)*l  ,017 

l6(jNfi.r,.Fn,ii  Oc  TO  “70  - 

l7fR(),GT.l,0fc5)  CO  TO  “80 

I'TLCoi  

papP*,00l*PP 

(Sal 

CO  TO  10 

COnTIRUE  

DtLXEXaXS 

I^f CELX.GT.O.O)  Y«»Y 

co  TO  eoo  

COniTJRUE 

THp  EOCTOP  I.OI7  IS  HfcRt  BECOUSt  OF  TnE  SLIGHT  ERROR  HETUHNEO 

tZ  from  syZYgy,  - 

iFtT.LT.PrOTTfFp.PS.NZn  Y»-Y 

CC  TO  800  

continue 

nFlOGiI  _ , . . 

CO  TO  800 

continue  . . - - 

nF|  tCst 

IF tlHETi.CT.2.617)  CO  TO  800 

nFi  tCiQ 

lFfTHLT4.LT. 0.3)  CO  TO  530  . __  . . 

TmFS®vncdm»( 1 .♦,5/rr«*3)»SIn(ThET*/2, J/RR*OTEl 

THnsTHETA"THFS  • - 

iF  f ThO.LT.O.Ol)  CO  TO  530 

XaRR»5lN(TH0) 

CO  TO  600 

continue  

pOiNT  lIfS  On  EXTREHELT  small  STREAM  FUNCTION  ANC  IS  CLOSE 
TO  RISE  AXIS--  approximate  X CCORD I N A TE t C*LCUL A TE  T FRUM  X, 


312 

xormaI ,01 

310 

CO  TO  600 

310 

5ao 

continue 

310 

CU  TO  800 

515 

560 

continue 

315 

TOiFsTZ-tLAST-T 

326 

X*SCRT(P).RLA33 

326 

yRl»*HS{RL2-XaX) 

331 

Y*-SQRT(yRL) 

337 

jF  (T,LT,Rn0TT(Fp,PS,N2)) 

350 

YaY*VNURM»TOlF 

353 

CO  .TO  800 

CIPHER. 168 
CIPHER. 16R 
CiPHf  R.1 70 
CIPHER. 171 
CIPHER. 172 
CIPHER. 173 
CIPHER. 170 

- cipher. 175 
CIPHER. 176 
CIPHER. 1 77 
CiPHtR.l 78 
CIPHER. 170 
CIPHER, 180 

--  cipher. 181 

CIphFR.182 
cipher. 183 
ClPHEP.180 
CIPHER. 185 

cipher. 186 

— CIPHER. 187 
CIPHER, 188 
CIPHER. 189 

In  CIPHER.190 
CIPHER. 191 
cipher. 192 
--  cipher. 193 
CIPHER. 190 
CIPHER. 195 
CIPHER. 196 
CIPHER. 197 
CIPHER. 198 
. CIPHER.  199 
CIPHER, 200 
CIPHER. 201 
CIPHER. 202 
cipher, 203 
CIPHER. 200 
CIPHER, 2U5 

cipher. 206 
CIPHER. 207 
CIPHER, 208 
CIPHER. 209 

cipher. 2l0 
CIPHER. 211 

cipher. 212 

CIPHER. 213 
CIPHER. 210 
cipher. 2l5 
CIPHER, 216 
ClPHtF.2l7 
CIPHER. 218 
_ CIPHER. 219 
CIPHER. 220 
CIPHER. 221 
CIPHER. 222 
CIPHER. 223 
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1 

i 


CIPHfP 

60(1  continue  cipher. 

^^U  - tiaf* ClPHth.^iiS 

5S5  ONnf‘  = 4BS(l(a-?,*P)*. 006001  CIPHER. 22H 

361  Y"*HS(  ()r?/0NpM).»,66666»*?)  - ClP“fR.?a7 

365  V*.SURT(yI  cipher. ^^8 

373  I^f T,LT.RnOTT(FP,PS.N2))  Y»»Y  - CIPHER. 22R 

1I0«  800  C0^T1NUE  ClPHtR.230 

ana  — rEturn cipher. 2ii 

*05  ERo  cipher. 232 


Subprogram  i fncth 
0163a 


) 


SCnCk 


SUrSOUTInF  ScHCKtXOU.K.ALPMj.xSiV)  -SChCk,? 

C SChCX.-i 

C THjS  SUBROUT|Kt  calculates  TmE  SnDCK  UISPLaCEMENT  SChCk,<i 

c OP  POINT  xoo  froh  event  k. giving  The  Nt*  Position  xsav,  schck.s 

c SCHCN.fc 

C inputs  - - SChCk.t 

C xCu  alNiTlAL  position  SCmCK.n 

C K»FVtNT  NllHflEP SChCn.r 

C ALPhAsShoCK  SCAliNS  parameter  SCmCk.jo 

C . - SCmCk.i  1 

C output  SChCk,i2 

C - SChCn.iS 

c xs*vsfin*l  Position  after  displacement,  schck.jr 

C SCmCs.iS 

C time  ntPENOCNT  hpOEL  parameters  CEcTD.2 

CCmmON  /gEuTo/  nF  . INOXF(SO)»  RTKbO).  HLFISO),  mK50)»  OCKSfl).  CEnTO.j 

1 CcFtSPIf  h»*XF(50)»  MMlNF(bO)»  K1NCF(S0)»  T I LTF  ( bO  ) , GFP.TO  , a 

2 AGt(50),  NpTA,  INOXO(IOO),  UtAoLdOO),  MUR(lOn),  GFoTO,? 

J HP|RC100)»  RTbsClon),  RLiJSdOu),  mMSCIOO#  HNBS  ( t 00 ) , GEpTO , 6 

« GCRTAf  t 00)  .r.LPTA  ( 1 00)  , IF  (SO) , TCh»R(N0)  ,mhGICCS0»  r GFdTD.T 

S XF»t50)»TFB(S0),Z(-»{S0),RUI(b0)  GEOTO.a 

COMMON/BAnnAT/XiN(J,)0),XtV(J,l0),RtlM10),RLvl1o5*VXyZ(J#l0),  R ado a 1,2 
t VRZ(|0)»CVLC(S»10)  RAOOaT.J 

C -----  HA[)DAT,(1 

c eventx  Parameters  evemx.f 

C^^mVON  /eventx/  NX,  lOX,  Tb(SO),  HfiCSO),  CCH(SO),  CLHISOJ,  - EVfMx.l 

1 lOGAp(bO),  PHCblSn),  MSP(bO),  TEmH(5o),  vP I SL ( 50 ) , E ve n t x . a 

2 RDZERO(SO),  M^ZFRnCbu),  kUZLkg(50),  BXo(50),  EVENTX.S 

1 BVStSO),  BZe(50)»  LHyBCSOJ,XALPMA(50),  KALCh  EveMX.6 

C . SCmCk.iR 

C FInaL  position  of  event  R,  first  find  distance  BETmEEN  SChCr,2() 

c event  and  point., SChCn.pI 

C SChCr,22 

dimension  XOuf J),*SAv(3),C(I),D(i) SChCN,2J 

CAlL  SuevFC(xOU,XlN(i,K),C)  SChCn,25 

II  SRsXMAG(C)  SCHCK.^b 

IT  SR|  aMbSR/alPm*  SCmCn,?T 

20  DLAM  = 0.0«fc/(SRL*M*(2,t8RL*M)J SChCK.25 

23  SRNE"aSR.OL*M*ALPHA  SChCk,2<> 

25  C*l  L UMtvIC.OI  - SChCx.jO 

27  CtlI*SRNFM»D(i)  SCmCk.jI 

SI  C(?)*SRNem«0(2)  --  - SCmCK,j2 

32  Ct5>®SRNEN«0(3)  SChCn.jS 

3«  aB,!  -SChCk.ju 

36  C*lL  VEClIN(aB,XINCI,K),A8.C,XSAV)  SChCk.jS 

«*  PEtURN  _ SChCK.jT 

«T  ENo  SCHCF.JB 


Subprogram  i fngth 

OOtOl 

function  asstCnmEnTS 
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I 

t 


"OOIT 


6 

10 

IT 

?t 

21 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


function  POnTT(FP,P3,N2) 

THjs  Function  supppoct***  cflcul*tes  tmc  tipe  contoub  that 

CROSSES  ZERO  »T  t (ilvEN  STNEAM  FUNCTION  VALUE 

INpUT  - - - 

INPUT  defined  in  cipher 


output 

HOnTT=TlME  Of  contour  CROSSING  ZERO  AXIS  FOR  ClVfN  X - - 
dimension  PS(iS),TCRS(15) 

data  TCBs/ii,IRb,3,M.2.T7.l.71,l,OR*,R8R».B0'».,56R»t«U8».503 

1 .22  I, . 12?7f  .068,  ,01SB»  tOOOli/ ^ 

AlsFP 

A2rl  .-FP  . - 

tC=1CRS(n2)»« A2*lC«StN2.1)*»Al 

RCnTTsYC  - - - 

return 

...end — 


POOTT,? 
HOOTT.3 
HOnTT  , (I 
RQOTT.S 
RnoTT,6 
HOoTT.T 
ROOTT.S 
ROoTT.R 
ROOTT,  to 
ROOTT.ll 
ROOTT, l2 
ROnTT,  !« 
ROoTT.iS 
— hoot  I. t 6 
ROOTT, )7 
HO(3TT,i8 
W"mT  . iR 
HOQTT ,20 
ROOTT, 23 
. -ROOTT, 2“ 


SUBRROr.RAM  I eNCTH 
000^6  . 


I 
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DISC** 

SUBROUTINE  0IsCtP(*«Q>D,H,SI>0IST,SR21)  - - - 

C 

c T-is  Routine  eincs  tkc  closest  dijt»nce  to'Tme  line  o*  centers 
C OR  INITUL  Event  RuSITicn  *>.0  E1n»L  EVENT  ROsIIlCNt 
C USES  vector  *n*lTS1$  *n0  ROSCCE  Ll^RSRTt 
C 

dimension  Pf3),o<S),C(J),R(S)fM(J)»M(J)  ' 

DImEvSIEn  d(J) 

*P«l 

12  C*LL  SuevECtO.P.C)  _ _ . 

15  C*LL  SuevECCO.P.F)  ■ • ' ‘ 

21  C*LL  PRCJ(E»C,CJ 

2*  e*LL  Su'JvECtS.F.R) 

11  DIST»»m*G(B) 

«2  ■ e‘LL  VECl1s(*P,P,*P»R.HT 

ah  SRaiPXMRCic) 

50  81«XMtc(C)  ■'  ■ 

c 

e DETERMINE  TNE  SENSE  0*  G RITM  RESPECT  TO  C — - - - 

C 

SJ *><»COT(C,G) 

0*  $)>$IGn($1,*M) 

T1  RETURN  ■■ 

71  END 


DISC**. 2 

0ISC*P,S 

01SC*P.R 

CIEC*P,5 

OISCAP.S 

DiSCAP.T 

UlSCiP.B 

01SC*P,« 

DlSCAP.IO 

OlSCAP.I I 

DISCAP.12 

DiSCtP.l) 

DlSC*P,l« 

DlSCAP.tS 

D1SC*P,|6 

Ol$C*P.i7 

DiSCAP.Ie 

DlSCAP.t9 

DI6CAP,20 

DlSCAP,2i 

DISCAP.22 

OISCAP.2} 

DlSCAP.Pu 

01$CAP,2s 


Subprogram  length 
00125 
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r 


EvCNiO 


lUsaouTiNE  EvrMo<$auo,^<i,pivz,iicu) 

Tt-IS  $i/8»bCcs»i  OCtS  T«e  vECTCR  AOOITIqn  Cf  THE 
TO  give  t FiNiU  OISPl'CEhEHT  

INPUTS 


EVtHPD,* 
EVtNAO, J 
PUSlTlOH  P*YZ  fVt>.AO,(i 
EVI«.A|J,5 
ETEMAD.t 


SPCC(1,L)iPOSITtCH  ARHay  OP  POINT  AFTER  AFpECI  PPOM  EVENT  I 
nAIpsUhBEb  vectors 

OUTPUT  ■ 

»CU(J)bFJNAL  vector  POSITION 


l« 

IP 

I* 
IT 
21 
21 
ST 
• I 
«1 
tT 
It 


too 


Dl«f nJtCn  SPCU(1,NaI),PXy2( J),»0U(3),C(J) 
call  »"IT(.J,o,0,CJ 

rb|ii,Ai.| 

00  >00  lBt,NH|  _ _ _ 

LL»L*1 

ABl 

B»1  • 

CALL  VECSUH(a,SpCU(I.I),3<SP0U(I>|.L)>C} 
CCnT;nuC 

rSUBoI  ,*FLOAY(igMl  ) 

call  VECLINCA,C,eSUB,PXY2,*0U) 

return 

END 


EVENAO.T 

EVENAD.S 

EVENAD.R 

EVENAD.IO 

evenao.ii 

EVENAO,I2 

EVI>.A0,1S 

EVENAD.IR 

fVENAO.lS 

EVENAO.l* 

EVENAO.IT 

EVtNAp.la 

EVF.NAO.IR 

EVENAO.ZO 

evenao.ii 

EVENA0,22 
fVEAA0.2I 
tvt>  AO.lB 
EVENAD.2S 
EVENAD,26 
EVENAO.27 

evenad.is 


IUBPRCCRAH  length 
00077  
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SECTION  VI 
TYPICAL  OUTPUT 

Three  test  problems  were  run  to  give  some  typical  output  from  the 
flow  field  models.  Each  test  problem  consisted  of  finding  a position  at 
burst  times  of  a given  particle  at  the  calculation  times.  The  first  two 
problems  consisted  of  a single  event  at  43  km  and  a spherical  Vortex  radius 
of  8.68  km.  The  last  test  problem  consisted  of  a five  burst  scenario  at 
relatively  small  yields  and  altitudes. 

The  output  gives  the  number  of  actual  events  that  have  occurred 
as  distinguished  from  hydromerged  events.  The  values  of  the  parameters 
are  given  for  both  the  initial  and  the  calculation  time.  The  initial 
parameters  are  labeled  with  the  "Type"  - INITIAL;  the  calculation  time 
parameters  have  a "Type"  - ACTUAL.  The  height,  colatitude  and  longitude 
are  given  along  with  the  indexes  of  the  fireballs  affecting  the  point.  The 
position  of  the  point  at  the  burst  time  is  then  given. 

For  example,  the  first  test  problem  has  a point  at  48  km  at 
2.7  sec;  at  0.0  sec  the  point  is  at  41  km.  The  drastic  change  is  altitude 
for  such  a short  time  span  is  a result  in  assuming  that  the  shock  arrives 
instantaneously  at  burst  time.  A correction  of  this  error  has  been  made 
and  will  appear  in  the  next  version.  The  shock  error  is  even  more  apparent 
for  the  given  point  at  70  km  at  2.78  sec  which  was  at  53  km  at  0.0  sec. 


wmmmm^rnmr 


OUTPUT  FOP  FLC^i  FIELD  Tf5T  ••  PPOBlEm  NUHBEP»  I 


event  no, 
\ 


time 

0,00 


NUMBEP  OF  ACTUAL  EVENTS  ■ I 
CCtWOI».ATES  AT  tUMST  TJMt  FOR 
MElCMTtuM)  CCLATITUCE  LUN'GITUDt 
U3,0000  Ml, 2066  230,0607 


I events 

RAUIuS(KM)  VtLOCITY(KM/S)  TVPE 
6,66003  OiOOOOOOOlNlTlAL 


coordinates  of  1 events  at  CaLCULATIUN  time  « 2,760 

event  no,  Tjme  heICNT(km)  cCLaTiTude  LUNGITUUE  HADIoSIKm)  VELOCITy(nm/S)  type 

1 2,76  03,1565  01,2066  236,0607  6,3o633  ,0362637  ACTUAL 

POSITION  OF  POINT  AT  CALCULATION  TIME  » 2,T6SEC 

MEICHT(Km)  COLATITUOt  LONGITUDE 

05,0000  01,2000  236,0607 


POINT  appears  inside  burst  number  I 


HEIGHT(Kn)  colatituoe  longitude 

31.0000  01,2000  236,0607 

number  of  events  affecting  PCInT  * I indexes  are  I 
POSITIONS  UF  puInT  at  event  TImES 

Event  time  mEIGhtifm)  cDlaTITUUE  longitude 
1 0,00  0O,30o0  01,2005  236,R8070 

MEIGNT(RM)  COLATITuDE  LQNCITuOE 

06.0000  01,2000  236,0607 

number  of  events  affecting  point  a I INDEXES  ARE  I 
positions  UF  point  at  event  times 

event  time  hEIOmTIOm)  colatituoe  longitude 
1 0,00  01,3261  01,2060  238,06070 

MEIGMT(Km)  colatituoe  longitude 

00,0000  01,2000  238,0607 


POINT  appears  inside  burst  NUMBER  I 


HElCMTtOM)  colatituoe  LONGITUDE 

70,0000  01,2000  23S,08o7 

number  of  events  affecting  PCInT  a I INDEXES  ARE  0 
POSITIONS  OF  point  at  event  TImES 

event  time  hEICmT(Km)  COLATITUOE  LONGITUDE 
t 0,00  53.3256  at.27Sl  23S.RS070 


COO»DIn»TES  Of  I EVr^TS  *T  CALCULitrCN  time  ■ 25.000 

Event  no.  time  meIC^-TIKM)  cCLATITUCE  LO^-I,ITUOE  HAi-IuSINM)  VELOCITy(KM/S>  type 

I 25,00  Sb.OSEe  ttt.EOSb  23d,96o7  Il.OPElT  .558tll«  aCtUaU 

POSITION  CE  POINT  *T  calculation  TIME  ■ 25,00SEC 

meicmt(km)  colatitude  longitude 

05. 0000  0I.24U9  230,9807 

number  of  events  affecting  point  • I INDEXES  ARE  I 

POSITIONS  OF  point  *t  event  times 

event  time  MEIGmTinm)  cOLaTITUDE  longitude 

I 0,00  44,2030  4I.248R  238,98070 

M£IGHT(Km)  COlATITUOE  LONGITUDE 

51.0000  4i,auu9  23S,Rao7 

POINT  appears  Inside  burst  number  t 

HEIGHT(Nm)  colatitude  LOnGITuDE 

48.0000  41,29119  236,9807 

POINT  appears  inside  burst  number  1 

MEIGnt(Km)  COLATITUOE  LONGITUDE 

44.0000  41,2449  236,9807 

Number  of  events  affecting  point  ■ i indexes  are  i 
POSITIONS  OF  point  at  EVENT  TIMES 

Event  time  meight(km)  colatitude  longitude 

1 0,00  43,3918  4I.244O  236.98070 

HEICmtikmJ  colatitude  longitude 

7o«0000  41.2440  236.46o7 

number  of  events  affecting  point  ■ 1 indexes  are  1 

POSITIONS  OF  point  aT  EVENT  TIMES 

Event  time  meigmt(km)  colatitude  longitude 

1 0,00  42.7415  41.2493  238.98070 
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coordinates  of  1 events  at  calculation  tine  ■ 50,000 

event  no,  TI^E  HEICnT(KM)  COLATITUCE  LONGITUDE  PAOIuS(K«)  VELOC I T Y (kM/3)  TVPE 

I 50,00  71,5TTT  ttl,2‘>88  23e.Nft07  I7,5«17<»  ,5715537  ACTUAL 

POSITION  OP  POINT  AT  CALCULATION  TIME  a bO.OOSEC 

MEICHTIMm)  COLATITUDE  LONGITUDE 
PS,0000  Vi, 2409  233,9607 

number  of  EVENTS  affecting  PCInT  • 1 INDEXES  ARE  I 

POSITIONS  OF  point  at  EVENT  TIMES 

event  time  MtlGMT(AMJ  colatitude  LONGITUDE 
t 0,00  4<i,7787  ai,2426  238.98070 

HEIOMT(Km)  colatitude  longitude 

51,0000  41,2<i«9  238,9807 


COQRDlNiTES  CF  } EVENTS  *T  C*LCUt*T10N  Tint  ■ 75,000 

CvtsT  NO,  time  HEIGHT(kM)  CCLATUcDt  tU^-CITuiiE  K*OIuS(ki)  VELOCITY  (kM/S)  TYPE 

1 75,00  83,6565  «l,2'»68  238,9807  22,69956  ,5«20e9l  ACTUAL 

POSITION  OF  POINT  AT  calculation  TINE  ■ 75,008eC 

HEICHTCKh)  COUATITUOt  LONGITUDE 

65.0000  61,2999  238,9So7 

nuhbeh  of  events  affecting  point  9 1 indexes  are  I 

POSITIONS  OF  point  at  event  times 

Event  TI>*E  HEICKT(Km)  cOLaTITUOE  longitude 

1 0,00  99,5750  91,2«13  238,98070 

NEIGHT(Km)  COUATITUDE  longitude 

51. 0000  91,2999  238,9807 

number  of  events  affecting  point  • 1 INDEXES  ARE  I 

POSITIONS  OF  POINT  AT  EVENT  TIMES 

Event  time  meigmt(km)  colatituoe  longitude 

I 0.00  99,9907  91.2511  238.98070 

HEIOHT(Km)  colatitude  LONGITUDE 

98.0000  91,2999  238,9807 

number  of  events  affecting  point  * 1 indexes  are  1 

POSITIONS  OF  POINT  AT  EVENT  TImES 

event  time  ME1GhT{Nm)  COlaTiTUOE  longitude 

1 0,00  97,0339  91.2953  238.98070 

meight(km)  colatitude  longitude 

99.0000  91.2999  238,98o7 

number  of  events  affecting  point  ■ 1 INDEXES  ARE  I 

POSITIONS  OF  POINT  at  EVENT  TImES 

event  time  meightimm)  colatituoe  longitude 

1 0.00  99.0560  91.2910  238,96070 

heiomt(kn)  colatituoe  longitude 

70.0000  91,2499  23S,9607 


POINT  appears  inside  burst  number  I 


OUTPUT  FOB  FuCx  FIELD  TEST 


PBO0LEM  NUPBEH* 


2 


i 

f 

[ 

E 

P 

I 

k 


NUMBEB  OF  actual  Events  p i 
CC'HOISATES  AT  buBST  Tl^E  FOB  J EVENTS 

EVENT  NO,  time  hElCiMtIKMj  cCLATlToDE  LONGITUDE  BADIlS(KM)  VELOCITY (aM/S)  TVPE 

I 0.00  UJ.OOOO  B1.2B8S  236,9807  6.66003  O.OOOOOOOINITIAI 


COOBoInATES  of  1 events  aT  calculation  time  • 25,000 

EVENT  NO.  time  heIghT(kM)  cCLATjTuct  LONGITUDE  UaDIuSIKm)  VELOC I T Y ( nM/S)  TYPE 

1 25,00  56.9528  91,2988  236,9807  11,69217  .5561119  ACTUAL 

POSITION  OF  POINT  AT  CALCULATION  TIME  ■ 25,00SEC 

HEicMTfKM)  colatitude  longitude 

03.0000  01,2986  236,9807 

NUMbFO  OF  events  affecting  POINT  • I indexes  ABE  1 
POSITIONS  OF  POINT  iT  EVENT  TIMES 

EVENT  TI-E  meiChT(Kmj  COLaTITUOE  LCnCITUDE 
1 0,00  38,3511  9l,2*>67  238,96096 

HflGMT(KM)  COLATITUOE  longitude 

01.0000  01,2099  236,9807 

NUN0EB  OF  EVENTS  AFFECTING  PoInT  ■ 1 INDEXES  ARE  I 

POSITIONS  OF  POINT  *T  EVENT  TIMES 

EVENT  time  hElGHT(AM)  COlaTtTUDE  LONGITUDE 

1 0.00  90.9063  91.2390  236.98070 

MEICMT(Km)  COLATITUDE  LONGITUDE 

39.0000  91,2999  236,9807 

number  of  events  affecting  POINT  » 1 indexes  are  1 

POSITIONS  OF  POINT  AT  EVENT  TIMES 

event  time  MEICHT(Km)  COLaTITUOE  LONGITUDE 

1 0.00  36.6152  91.2371  236.98070 

meight(kh)  colatitude  longitude 

35.0000  91,2999  236,9807 

number  of  events  affecting  point  ■ 1 indexes  are  I 

POSITIONS  OF  point  at  event  times 

event  time  m£ICmT{«m)  cOlaTITUDE  LCNCITUDE 

1 0,00  35,0228  91,2381  238,98070 

meight(km)  colatitude  longitude 

25.0000  91,2999  238,9807 

NUMBER  OF  EVENTS  AFFECTING  PcInT  • 1 INDEXES  ARE  0 

P08ITIONS  OF  POINT  aT  event  TIMES 

event  time  meigmt(km)  colatitude  longitude 

1 0.00  25.2216  91.2950  236.98070 
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cr>o«oiN*Tes  of  i events  *t  cauculation  tihe  ■ so.ooo 

event  no,  time  HEIChT(km)  cCLATITuDE  LOnCITUUE  HAOIuS(KM)  VEUOCITy(KM/S)  type 

I so, 00  71,5777  «l,2V«8  210,9807  17,1«179  ,5715517  ACTUAL 

POSITION  OF  POINT  AT  CALCULATION  TIME  ■ SO.OOSEC 

heigmt(km)  colatitude  longitude 

PS, 0000  91,2986  210,9807 

NUHbEB  op  events  affecting  point  » 1 indexes  are  1 

POSITIONS  OF  point  at  event  times 

event  time  MEIGmT(Km)  colatitude  longitude 
1 0,00  18,1162  al,2987  218,98099 

HEICHi(hh)  colatitude  longitude 

91.0000  91,2009  218,9607 

number  of  events  affecting  PCInT  a 1 INDEXES  ARE  1 
POSITIONS  OF  point  at  event  TIMES 

event  time  MEIGmT(km)  colatitude  longitude 

1 0,00  40,6565  91.21PP  2‘i6,'*6070 

MEICHT(Km)  colatitude  longitude 

19.0000  91,2009  218,9807 

number  of  events  affecting  point  9 1 indexes  are  1 

POSITIONS  OF  point  at  EVENT  TImES 

event  time  HEICmT(Km)  colatitude  longitude 

1 0,00  18,5610  01,2155  218,98070 

HEICHT(Mm)  colatitude  longitude 

15.0000  91,2009  210,9807 

number  of  events  affecting  point  p 1 INDEXES  ARE  1 
POSITIONS  OF  point  aT  EVENT  TIMES 

event  time  HEICMT(Km)  colatitude  longitude 

1 0,00  19,7814  01,2175  218,58070 

HEIGMT(Km)  colatitude  longitude 

25.0000  91,2049  210,9807 

number  of  events  affecting  point  b 1 INDEXES  ARE  0 
POSITIONS  OF  POINT  at  EVENT  TIMES 

event  time  meight(km)  colatitude  longitude 

I 0.00  25.2218  01.2950  238,98070 
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C0ni»0lN*TE8  OF  1 FVpNTS  AT  CALCULATION  TIME  ■ 75,000 

event  no,  time  HLIGMT(kM)  cOLATiTUOE  LJNGITUUfc  MADIuSikm)  VtLOC I T Y ( KM/S)  TVPE 

1 75,00  ei,b56i  <ll,f9a8  iijs.oeo?  ,>2,b94l5t>  , 5^20841  ACTUAL 

POSITION  OF  POINT  AT  CALCULATION  TIME  m 75,00SEC 

HEXCHTIKm)  COLaTITUOE  LONGITUDE 

43.0000  4t,a9SS  238,9807 

NOPBER  of  events  affecting  point  ■ I INDEXES  ARE  t 
POSITIONS  OF  POINT  at  event  TIMES 

event  time  MEIGHtfKM)  COLATITUOE  LONGITUDE 

1 0,00  38,3387  41,2487  236,46095 

HEICHT(Km)  COLATITUOE  LONGITUDE 

41.0000  41.2990  238,4607 

NUMBER  OF  EVENTS  AFFECTING  POINT  » 1 INDEXES  ARE  I 

POSITIONS  OF  PUInT  at  event  TIMES 

event  time  MEIGHTIKm)  COLATITUDE  LONGITUDE 

1 0,00  90,1493  91.2370  236.46070 

height(km)  colatitude  longitude 

39.0000  41,2499  238,4607 

number  of  events  affecting  point  a I INDEXES  ARE  1 
POSITIONS  OF  point  at  EVENT  TIMES 

event  time  MEICHT(XM)  COLATITUOE  LONGITUDE 

1 0,00  36,2460  41,2356  236,46070 

MEIGHT(Km)  COLATITUDE  lonsitude 

35.0000  41,2944  238,46oT 

number  of  EVENTS  affecting  POInT  ■ 1 INDEXES  ARE  1 

POSITIONS  OF  POINT  at  event  TIMES 

event  time  MEIGMTIKmI  COLATITUOE  LONGITUDE 

1 0.00  34,9234  91.2374  238,46070 

MEIGMTIKM)  COLATITUOE  LONGITUDE 

25i000O  41,2944  238,4607 

number  of  events  affecting  point  ■ I indexes  are  0 

POSITIONS  OF  point  at  EVENT  TIMES 


event  time  HEICMT(Xm)  COLATITUOE  LONGITUDE 
1 0*00  25,2216  41.2450  236,46070 
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OUTPUT  POP  PLC"  FIELD  TfST 


problem  NlirBfcRi 


) 


event  no. 

time 

nUmpeR  op  actual  events  » 5 

COORDInatfs  at  bijRST  TImF  FOR 
MEIGhTIkM)  CCLATITUOE  LONGITUDE 

5 events 
HADIuSCKfl) 

velocity(pm/s)  type 

1 

0.00 

5.0000 

Rl,29aB 

2J9,009J 

,2bT05 

O.OOOOOOOlNiTlAL 

2 

0.00 

R.lSno 

ttl,298P 

236,980T 

,2b705 

O.OOOOOOOInITUL 

s 

.20 

R.5000 

Rl ,2968 

238,9807 

,26705 

O.oOoOOOOlNlTlAL 

<1 

12.00 

12.0000 

01.2988 

239,0093 

.26705 

O.OOOOOOOInItIaL 

5 

20.00 

11.5000 

41,2988 

239,0093 

.26705 

O.OOOOOOOlNiTlAL 

COORDIn 

ATES  OF  5 

events  aT  CALCULATION  TIME 

a 22,000 

event  no. 

time 

Mtl&MT (K>0 

CCLATI’L'OE 

longitude  RAOIuS(KM)  VELOCITyCKM/S) 

TYPE 

1 

22.00 

5,9868 

41,2986 

239,0093 

,55334 

,0448540 

ACTUAL 

2 

22.00 

10,2207 

41,2988 

238,9807 

,52425 

,0486689 

ACTUAL 

3 

22.00 

10,6602 

41,2988 

238,9806 

.53710 

.0532197 

ACTUAL 

4 

22.00 

12.7610 

41,2988 

239,0093 

,56068 

,0761038 

ACTUAL 

5 

22.00 

11,50B3 

41,2988 

239,0093 

,80005 

,0041602 

ACTUAL 

POSITION  OF  POINT  at  CALCULATION  TIME  c 22,005eC 

HEICHT(Kh)  CCLATITUDE  longitude 
R.OOOt)  «l.2<»be  238,<»8b<» 

number  of  FVEnTs  aFPECTIng  point  a 1 INDEXES  ARE  2 
POSITIONS  OP  POINT  aT  event  TIMES 


event 

time 

meiGwT(Amj 

COLaTITUOE 

LONGITUDE 

1 

0.00 

8,9472 

41.2988 

238,98617 

2 

0,00 

8,9509 

41.2988 

258,98615 

3 

.20 

8,9539 

41.2986 

238,98650 

4 

12,00 

8,9343 

41.2988 

238,98650 

5 

20,00 

8,9745 

41,2988 

238,98679 

heigmt(fm)  colatitude  longitude 

9.1500 

41,2968  238,9864 

Number 

OF  EVENTS 

affecting  PCInT  a 

2 indexes  are  2 

POSITIONS  OP  POINT  aT  event 

TIMES 

Event 

time 

mEIGmT(am) 

COLATITUDE 

LONGITUDE 

1 

0.00 

9,0902 

41.2988 

236.98613 

2 

0.00 

9,0938 

41.2986 

238,98611 

3 

.20 

9.0276 

41,2968 

238,98677 

4 

12,00 

9.0163 

41,2988 

238,98656 

5 

20.00 

9,0570 

41,2988 

238,98662 

neicht(nm)  colatitude  longitude 

9,9o00 

41.2988  238,9864 

Number 

OF  EVENTS 

affecting  PCInT  a 

2 indexes  are  2 

POSITIONS  OF  POINT  AT  event 

TIMES 

Event 

time 

MEICMT(aM) 

COLaTITUDE 

LONGITUDE 

1 

0.00 

9,9259 

41,2988 

238,98459 

2 

0,00 

9,9289 

01.2986 

238,98457 

3 

a20 

9,8821 

01.2968 

238.98479 

4 

12,00 

9,9681 

41.2988 

238,98715 

s 

20,00 

9,6951 

41,2966 

238,98776 
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I 

j 


i 


HElGHTtKM)  C(jL*T1TUDE 
10,0800  «1.2900 


LONGITUDE 

^30,9e6« 


POINT  APPEARS  INSIDE  BURST  NUMBER  } 


MEICHT(Rm) 

11,2000 


COL*TITUDE 
«l  2980 


LONGITUDE 

210,9864 


number  or  EVENTS  AFFECTING  PCInT  ■ 1 INDEXES  ARE 

POSITIONS  OF  point  aT  EVENT  TIMES 


event 

time 

height  cKm) 

COLATITUOE 

LONGITUDE 

1 

0,00 

10.9509 

01.2968 

238,90007 

2 

0.00 

10.9615 

01.?986 

238.98006 

3 

.20 

10.9726 

01,2988 

238.98009 

0 

12,00 

11.0068 

01.2968 

238,90505 

5 

20.00 

11.1508 

01,2966 

238,90800 

COPPOINATES  OF 

5 events  aT 

CALCULATION 

30,000 


NO. 

time 

MEIGMT (KM) 

CCLATlTUOt 

LONGITUDE 

HADIUS(KM} 

VELOCITy(KM/S) 

type 

1 

30.00 

6,3805 

01,2960 

239,0091 

,66313 

.0433500 

ACTUAL 

2 

30,00 

10.5020 

01.2980 

218,9600 

,65083 

,0«6O166 

ACTUAL 

3 

30.00 

10.9900 

01,2988 

230,9003 

,60061 

,0500237 

ACTUAL 

0 

30.00 

11.3330 

41,2988 

239,0091 

.67395 

,07o0575 

ACTUAL 

5 

30,00 

12,2092 

41,2968 

239,0093 

1,06126 

,07«9i70 

ACTUAL 

POSITION  OF  POINT  AT  CALCULATION  TIME  s 10,008EC 

MElGhT(KM)  COLATITUOE  LONGITUDE 
9,0000  41,2980  210,9064 


number  of  events  affecting  point  » 1 indexes  are 

POSITIONS  OF  POINT  at  event  TIMES 


event  time 
1 0.00 
2 0,00 
S .20 

4 12,00 

5 20f00 


MEIGmT(KM) 

8,9360 

8.9198 

0,9U2T 

0.925T 

e.96|a 


COLATITUOE 

01,2986 

01.2966 

41.2908 

01.2988 

01.2968 


LONGITUDE 

230,90620 

218.96617 

238.98652 

238.98653 
236,98601 


MEICMT(Km)  COLATITUOE  LONGITUDE 
9.1500  41,2906  230,9064 


number  of  EVENTS  AFFECTING  PCInT  » 2 INDEXES  ARE  2 S 

POSITIONS  OF  PUInT  at  EVENT  TIMES 


event 

time 

HE1GHT(K«.) 

colatttude 

LONGITUDE 

1 

0,00 

9.0802 

01,2900 

236,90616 

2 

0.00 

9,0830 

01.2908 

238.90610 

3 

.20 

9.0038 

01,2960 

230,90602 

4 

12.00 

0,9956 

01.2960 

230,90662 

6 

20.00 

9.0300 

01.2900 

230,90607 
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HEICHT(km) 

COLATITUOE 

longitude 

9t9o00 

«l,29eA 

230.9669 

number  of  events  affecting  POINT 

■ 2 indexes 

ARE  2 

s 

POSITIONS  OP  point  *T  EVfcNT  TIHtS 


event 

TIME 

HEIGHT(KM) 

COLaTITUDE 

longitude 

1 

o.nn 

10.0922 

UI.2968 

?38. 98593 

2 

0.00 

10.0as2 

Ul.2986 

238.98592 

3 

.20 

10.0131 

fll.29«8 

238.98637 

4 

12.00 

10.0539 

01.2988 

238.98801 

5 

20.00 

10.0126 

01.?9b6 

238,98916 

HEICHTtPM)  COUATITUOE  LONGITUDE 


lOtBOOO 

91,2988 

238.9869 

POINT 

APPEARS 

INSIDE 

BURST 

number 

2 

POINT 

APPEARS 

INSIDE 

burst 

number 

I 

HEIGHTTPm)  COLiTITOOE 
lliZOOO  «1.2«e8 


LONGITUDE 

23e,<)66a 


POINT  appears  inside  burst  NUHBER  S 
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I 


event  no. 

coordinates  of  fe 

time  MtltNT(KM) 

events  at  calculation  time  s 40,000 

CCLATIfUOE  LUnUITUDF  RAOIuS(Kh)  V LLOC 1 T T Um/S) 

type 

1 

40,00 

6.6980 

41 ,29AS 

2J9,0095 

.78772 

,0424488 

ACTUAL 

2 

40,00 

JP,9J61 

41 ,29flB 

218,9804 

.75159 

.0446572 

ACTUAL 

3 

40,00 

11.1972 

41 ,29SR 

218,9802 

,76687 

,0476776 

ACTUAL 

a 

40,00 

13.8771 

41,2988 

23'>,0093 

,81183 

,0670375 

merged 

5 

40,00 

13.1745 

41 ,2988 

239.00R3 

1,18871 

,0837271 

merged 

6 

40,00 

12.7447 

41,2988 

239,0093 

, 10076 

o.ooooooo 

merged 

POSITION  OF  POINT  *T  CALCULATION  TIME  z <I0,00SEC 

HEIGHT(Km)  COLATITUOE  LONGITUDE 
9i0ooo  Mi.29ee  2}<i,9eb4 


number  of  EVENTS  affecting  PCInT  « I INDEXES  ARE 
POSITIONS  OF  POINT  at  event  TIMES 


event 

time 

MEIOHT(KM) 

COLaTITUOE 

LONGITUDE 

1 

0.00 

8. 9217 

41.2986 

218,98621 

2 

0,00 

8,9255 

41.2988 

216,98618 

3 

• 20 

8,9261 

41.2968 

216,98652 

6 

40.00 

8,9083 

41,2988 

218,98654 

5 

40,00 

8,9446 

41,2968 

218,98681 

heicht(km)  culatitude  longitude 

9.1s00 

91.2988  218.9869 

number 

OF  EVENTS 

AEEECTINC  point  b 

2 indexes  are  2 

POSI 

TIONS  OF  point  at  event 

TIMES 

event 

TIME 

MEIChtCXm) 

COLaTITUOE 

LONGITUDE 

1 

0,00 

9.0677 

41.2988 

218,98617 

2 

0,00 

9,0714 

41,2988 

218,98615 

3 

.20 

8,9734 

41,2988 

218,98684 

4 

40.00 

8,9614 

41,2986 

218,98661 

S 

40,00 

8,9999 

41.2968 

238,96688 

M£1GMT(Nm)  colatitude  longitude 

9.9000 

91.2988  238.9869 

Number 

‘ OF  EVENTS 

aFFFCTIng  POINT  ■ 

2 INDEXES  ARE  2 

POSITIONS  OF  point  at  event 

TIMES 

Event 

TIME 

meICmT(Xm) 

COLATITUOE 

LONGITUDE 

1 

0,00 

9,6l6l 

41,2968 

238,98655 

2 

0.00 

9,8194 

41.2988 

238,98654 

3 

.20 

9,8084 

41,2968 

238,98696 

4 

40.00 

9,7927 

41.2968 

238,98927 

5 

40.00 

9,7619 

41,2988 

236,98903 
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